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1. INTRODUCTION 

Supersymmetry [Biai3llia[ai2l[8lll[l0l[niIiaiI3IIllI3[ia[IIlIlH]isoneofthc most appealing concepts of 
physics beyond the Standard Model (SM). Some of the main motivations for studying supersymmetry are 

• Supersymmetry is an extension of the Poincare algebra which relates fermions to bosons. 

• As a local symmetry, supersymmetry is naturally connected to gravity. 

• Supersymmetry stabilizes the hierarchy [T51 [TOl EO] between the electroweak and the GUT or Planck scale. 
Thus, the quadratical divergences occurring in the Higgs-boson mass loop corrections are systematically canceled 
and fine tuning between the bare Higgs-boson mass and the loop contributions avoided. 

• Supersymmetry, realized below the TeV scale, unifies the SU{3)c x SU{2)l x U{1)y couplings at the GUT 
scale. 

• Supersymmetry provides a cold dark matter candidate, supposed matter-parity is conserved. 

A minimal supersymmetric extension of the Standard Model (MSSM) was proposed already some time ago. For a 
review of the MSSM see for instance Martin [5T] and references therein. In the MSSM the minimal particle content 
is added to the SM in order to arrive at a supersymmetric model. The MSSM extents the SM by an additional 
Higgs doublet, necessary to give masses to up- and down- type fermions and in order to keep the theory anomaly 
free. Every field is promoted to a superfield, pairing fermionic and bosonic degrees of freedom. However, none of the 
additional predicted particles has been observed up to now. Thus the question arises, why should we study further 
extensions of the minimal supersymmetric model? The reason is, that the MSSM does certainly not parameterize 
any supersymmetric extension of the SM. This is in particular due to the fixed particle content we encounter in the 
MSSM. For instance, although the Higgs sector in the MSSM consists of two doublets in contrast to one in the SM, 
the Higgs sector turns out to be highly restricted. Thus, at tree-level, the lightest CP-even Higgs-boson is predicted 
to be lighter than the Z-boson. Large quantum corrections are necessary in order to comply with the experimental 
lower bounds from LEP on the CP-even Higgs-boson mass, which in turn require a very large scalar-top mass. That 
is, some new kind of fine-tuning is to be introduced in the MSSM. 

Let us collect some reasons, making it worthwhile to study extensions of the minimal supersymmetric Standard 
Model: 

• In the MSSM we encounter the so-called /i-term in the superpotential, where /i is a dimensionful parameter. 
This parameter has to be adjusted by hand to a value at the electroweak scale, before spontaneous symmetry 
breaking occurs [22 . This is seen to be a problem of the model. It is desirable to look for further extension, 
which do not have this /^-problem. In singlet extensions, for instance, an effective /i-term may be generated 
dynamically. 

• As mentioned before, the Higgs sector is highly restricted in the MSSM. The lower bounds on the Higgs-boson 
masses from LEP measurements require large quantum corrections accompanied by a large stop mass in this 
model. An extended Higgs sector may relax this restrictions and thus circumvent the lower experimental bounds. 

• The MSSM Higgs-boson sector is CP-conserving at tree level. Extending the Higgs sector in an appropriate 
way, CP violating phases arise. Sufficient CP- violation would meet one of the necessary Sakharov criteria |23j 
in order to generate the baryon-antibaryon asymmetry in our Universe. 

• The baryon-antibaryon asymmetry may be generated by strong electroweak phase transitions of first order |24l 
l25l |26| [27] . The required cubic terms in the effective potential arise in the SM and the MSSM only via generically 
small radiative corrections. An explicit cubic term is possible in further extensions of the MSSM. 

Here we will review the next-to-minimal supersymmetric extension of the Standard Model (NMSSM) [25", [501 
[3Il[31l3Sll31[351l3S[371[MllSn], which has the capability to solve the mentioned hmitations of the MSSM. In the 
NMSSM an additional gauge singlet is introduced which generates the /x-term dynamically, that is an effective /i-term 
arises spontaneously and the adjustment by hand drops out. This is surely the main motivation for the NMSSM and 
may justify the price to pay, that is the introduction of an additional gauge-singlet superfield. The particle content 
in the bosonic part of the singlet results in two additional Higgs-bosons whereas in the fermionic part we have one 
additional neutralino, called singlino. Altogether we have seven Higgs bosons and five neutralinos in the NMSSM, 
compared to five Higgs bosons and four neutralinos in the MSSM. The Higgs-boson sector of the NMSSM is no longer 
CP-conserving at tree level, merely CP-conservation only arises if the parameters of the Higgs-boson sector are chosen 
in an appropriate way. Nevertheless, in most studies in the literature the special case of a CP-conserving Higgs sector 
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is considered, in order to simplify matters. In case of a CP-conserving Higgs sector we encounter altogether three 
CP-even Higgs-bosons, two CP-odd ones and in addition two charged Higgs bosons in the next-to-minimal model. As 
we will see, the Higgs-boson sector is in deed much less restricted and the lower mass bound prediction of a CP-even 
Higgs boson in the MSSM is generally shifted substantially. The Higgs-boson phenomenology can in general be 
very different from what to expect in the MSSM; in addition to supplement the total number of Higgs-bosons. For 
instance, in the NMSSM the possibility arises that a CP-even Higgs boson decays into two very light CP-odd ones 
which would have escaped detection at LEP and may even be difficult to detect at the LHC. Also the Higgs potential 
is enriched, leading to interesting consequences. Let us also address the last aforementioned advantage of the NMSSM 
over the MSSM. As we will see, the trilinear A-parameter soft supersymmetry breaking terms, corresponding to the 
superpotential, may account for the desired strong first order electroweak phase transition without large fine-tuning. 

The additional neutralino, that is the singlino, in general mixes with the other four neutralinos. Also in the 
neutralino sector there may be a substantial change of phenomenology compared to the minimal supersymmetric 
model; in addition to supplement the total number of neutralinos by a fifth neutralino. This is due to the fact 
that the singlino is introduced as a gauge singlet. Only through mixing with the other neutralinos this singlino 
has couplings to the non-Higgs particles. This opens the intriguing possibility to have a singlino-like neutralino 
which moreover may become the lightest supersymmetric (partner-)particle (LSP). Relying on matter-parity, every 
supersymmetric partner particle will in this case eventually decay into this singlino-like LSP. In particular the 
next-to-lightest supersymmetric (partner-)particle (NLSP) would, due to the small couplings, decay very slowly 
into the LSP. These large decay length' could be revealed by signatures of displaced vertices in the detector. But 
of course, in case the singlino-like neutralino is not the LSP, it would be omitted or at least be suppressed in 
cascade decays. Generally, in discussing the NMSSM phenomenology we draw the attention to differences which 
may occur in contrast to the MSSM. For the MSSM phenomenology we refer to the extensive literature to this subject. 

We start in Sect. |2] with the NMSSM superpotential. We review the main motivations for the modifications 
compared to the MSSM and also discuss the drawback which arises in context with the Pcccci-Quinn symmetry, which 
is promoted to a discrete Ka-symmetry by the introduction of the singlet cubic sclfcoupling in the superpotential. 
We recall the arguments in order to circumvent the occurrence of dangerous domain walls which emerge in context 
of spontaneously broken discrete symmetries. 

In Sect. |3] we present the mass matrices and the parameters to describe the complete Higgs-boson sector at 
tree-level. In this we treat the most general Higgs sector with the possibility of CP-violation. The special case of a 
CP-conserving Higgs sector may be easily inferred. We briefly mention the one-loop effective potential involving the 
ficld-dcpcndcnt Higgs-boson masses. Then we discuss the Higgs-boson phenomenology, where we stress the prospects 
at the LHC. In particular we review the discussion of a "no-lose" theorem, that is the interesting question, whether it 
can be guaranteed that at least one Higgs boson of the NMSSM will be detected at the LHC - supposed the NMSSM 
is realized in Nature. 

In Sect. [4] follows a recap of the neutralino sector, where we present the mixing matrix as well as neutralino 
phenomenology, which is of special interest in case there is a singlino-like LSP. 

In Sect. |5] we consider parameter constraints in the NMSSM, beginning with theoretical constraints. For instance 
a theoretical constraint comes from the requirement to have a potential with a global minimum not breaking electric 
or color charge. A further constraint restricts coupling parameters by forcing them to be pcrturbative up to the 
GUT scale. Of course the latter constraint relies on the fact that the model is valid in the large range up to the 
GUT scale. We will also discuss briefly an approach to quantify fine-tuning in general extensions of the SM. On 
the experimental side we start with considering constraints coming from colliders. Precision measurements of the 
.Z-boson width are discussed, followed by exclusion limits from searches for neutralino and chargino pair production, 
as well as Higgs-boson production. The possible contribution of the NMSSM to the anomalous magnetic moment of 
the muon is discussed. The constraints from the — > 57 decay are presented, which is a loop induced process and 
thus highly sensitive to possible new particles emerging in the loops. Turning to cosmological constraints we start 
with the very strong constraint origination from the new WMAP result on the relic abundances of the LSP, supposed 
this is the candidate for the observed cold dark matter. The direct dark matter detection experiments are recalled 
briefly. Then we come to the prospects of the NMSSM in terms of strong first order electroweak phase transitions in 
order to explain the baryon-antibaryon asymmetry. Finally, we consider some of the newer parameter scans, which 
provide some interesting insight for rather large ranges of parameter space. 
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In section [6] eventually we inspect the Higgs potential with respect to stationary solutions. A method is introduced 
employing a Groebner basis approach in the framework of gauge invariant functions, which allows to certainly 
determine the global minimum in an algebraic way and reveals a quite surprising stationarity structure. However, 
this method is restricted to the tree-level potential. 

We close with a rather extended appendix. We mention some of the frequently used computer tools, mainly used 
in parameter scans, followed by a derivation of the essential new Feynman rules in the NMSSM compared to the 
MSSM. Then we present the formalism of gauge-invariant functions as well as the Buchberger algorithm which allows 
to compute a Groebner basis, employed in the algebraic determination of the global minimum. 

In advance we apologize for any missing references in the bibliography. But in view of the vast amount of publications 
with view on the NMSSM, in particular in the recent years, this seems to be hardly avoidable. 
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TABLE L Particle content of the NMSSM in terms of decomposed supermultiplets together with the corresponding multi- 
plicity respectively hypercharge with respect to SU{'i)c x SU{2)l x U{1)y- Compared to the particle content of the minimal 
supersymmetric extension the new ingredient in the next-to-minimal extension is the Higgs supermultiplet S. Only the first 
family is given; the other families are introduced analogously. The here used convention for the hypercharge is Q = 1^^ + Yw 
with Q the electromagnetic charge in units of the positron charge, I^r the weak isospin 3-component, and Yw the hypercharge. 



2. THE NMSSM SUPERPOTENTIAL 



In the NMSSM the superpotential is 

W = %„(QTei/„) - dyaiOTeHa) - eyeCL^eHa) + \S{H^eHd) + ^^5^ (2.1) 

with dimensionless couplings yu, yd, ye, A and k. Note that the superpotential has cubic mass dimension. The 
supermultiplets are denoted with a hat and are given in Tab. |l] together with its bosonic and fermionic particle 
content and its multiplicity respectively hypercharge with respect to SU{3)c x SU{2)l x U{1)y- Compared to the 
MSSM particle content there is only one new supermultiplet ingredient, namely the gauge singlet S, that is one 
complex spin-0 singlet S and one spin-1/2 singlino S. The weak isospin indices are not written explicitly and e is 



defined below in (3.3 1. Note that the convention for the signs of the A and k terms in the superpotential differ in the 
literature. 



2.1. Discussion of the superpotential 



We start from the MSSM potential in terms of scalar fields, 

W^MSSM = uy^iQ^eH^) - d*^yd{Q^eHd) - e*j,y,{L^ eH^) + KH^eHd) . (2.2) 

This is the minimal superpotential, where we have instead of the dimensionless terms AS* (^H^eHd) + ^f^S^ a so-called 
/i-term, with jj, a dimensionful parameter. In any case a term proportional to H^eHd is required to give masses to 
up- and down-type fermions. This follows immediately from the F-terms derived from the superpotential; as given 



in (Dill. On the other hand we know from experiment that the Higgs vacuum-expectation value is ?; ~ 246 GeV, 
that is of the order of the electroweak scale. With view on the F-terms in the MSSM potential we get potential terms 
(l/ip -I- rn^^) l-f^uP + + '^Ifd) l-^d (This can be seen in the following discussion of the Higgs potential and can 



be read off from (3.4 1 and (3.7 1 if we replace jApISp by |/J,p). With the vacuum-expectation- values (Hd) = Vd/V^, 
= w„/\/2 and v"^ = + we see that without fine-tuning between the mass- and /^-parameters, the parameter 
yU can not be much larger than v. This is, that we have to adjust the /z-parameter by hand to the electroweak scale 
which in the SM arises from spontaneous symmetry breaking. This is the so-called ^-problem of the MSSM. 



1. Peccei-Qmnn symmetry 

The main idea is to circumvent this dimensionful parameter /i by introducing a new scalar field S which gets a 
vacuum expectation value to generate this parameter effectively, that is the /i-term is replaced by AS* {H^eHd) and 
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the required /x-term at the electroweak scale arises from spontaneous symmetry breaking with = X{S). In this 
way we have achieved a quasi-universal mechanism of spontaneous symmetry breaking. This is considered to be the 
main motivation of the NMSSM. For some other solutions to circumvent the /x-problem we refer the reader to Refs. 
[221 SOI mi BU 1131 HH SB SSI HZ] • A comparison of the Higgs sector of several non- minimal supersymmetric extensions 
can be found in Ref. |48j Now, since the bilinear /i-term in the superpotential is replaced by a trilinear term, we can see 
that the superpotential is invariant under a global phase transformation, the so-called Peccei-Quinn (PQ) symmetry 
U{l)pQ [inilSir. If we assign PQ charges qpQ to the chiral supermultiplets, denoted by 0,^ with i running over all 
chiral supermultiplets, in the following way [51 , 
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(2.3) 



we see immediately that the superpotential is invariant under the global U{1)pq transformations 

4)i -> exp{iqpQ9)4>i , (2.4) 

with an arbitrary phase 6. Note that since we assign PQ-charges to supermultiplets also the derived Lagrangian stays 
invariant under U{1)pq. In the minimal supersymmetric extension, MSSM, the PQ symmetry is explicitely broken 
by the ^-term. In the NMSSM this PQ-symmetry is spontaneously broken once the Higgs bosons acquire vacuum- 
expectation-values. This spontaneous broken continuous symmetry gives, following the Goldstone theorem, inevitable 
a massless so-called Peccei-Quinn mode. Such axion has not been found and exclusion limits were derived [52] . 
Experimentally there remains only a rather small corresponding A-parameter occurring in the superpotential term 
XS{H^ eHd) in the range 10""^ < |A| < 10""'^°. In order to get an effective /i term of the electroweak order, the 
vacuum-expectation-value (VEV) of S would have to be much larger compared to the electroweak scale - again a 
fine-tuning would be required. 

In the NMSSM thus an additional term is imposed in the superpotential in oder to violate the Peccei-Quinn 
symmetry by means of a cubic selfcoupling term, ^kS^, with k the PQ-symmetry breaking parameter. In this way 
all parameters of the superpotential are still dimensionless. Alternatively, one could impose a linear or quadratic 
term in order to break the Peccei-Quinn symmetry but then again a dimensionful parameter is necessary. Note 
that the superpotential has to have mass dimension three. Of course any term higher than trilinear in the fields in 
the superpotential is forbidden by the requirement of renormalizability. Finally we arrive at the proposed NMSSM 
superpotential, written in its scalar form 

W = uyuiQ'^eHu) - d*jiydiQ^eHd) - e*iiy,{L^ eHa) + XSiH^eHd) + ^kS"" . (2.5) 



2. Discrete TLz symmetry 



However, there is some debate, whether the NMSSM superpotential (2.5 1 is viable. The main concern arises since the 
continuous Peccei-Quinn symmetry is broken explicitly by the \k,S^ term, but there remains a discrete symmetry 
of the superpotential [55] : 

^e'^^;, (2.6) 

where the index i runs over all superfields occurring in the superpotential. This discrete symmetry is spontaneously 
broken by the VEV of the complex scalar field S. Inevitably, this spontaneous breaking of a discrete symmetry leads 
to the domain wall problem, a special kind of a topological defect. During electroweak phase transition of the early 
Universe, this spontaneously broken discrete symmetry causes a dramatic change of the evolution of the Universe and 
would spoil the observed cosmic microwave background radiation. 

A study of the generation of domain walls in the NMSSM potential with numerical methods can be found in 
Ref. [54 . Following Refs. [53l l55] we want to sketch the mechanism of the formation of domain walls in a model 
with spontaneous breaking of a discrete symmetry in the simple Goldstone model of a complex scalar field with 
Lagrangian 

C=l{d,<P^){d'^<j>)^^{<t>'~rj')\ (2.7) 

where 77 is a constant. The potential of this Lagrangian is the double well potential. The minima of the potential are 
at 



(2.8) 
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that is we have in this case a discrete TLi symmetry of the Lagrangian which is spontaneously broken by the vacuum. 
Now, the domain wall is nothing but the intermediate region between these degenerate minima (^±. Let us assume, 
without loss of generality, that this domain wall is situated in the y-z-direction. From the Goldstone-model Lagrangian 
we determine the Euler-Lagrange equation + {<^^ — rf^ = 0. The solution of this equation is 

= r,tanh(|-), (2.9) 

with thickness of the wall A = (A/2)^^/^ry^^. This solution describes the field in the intermediate domain wall region 
between the two vacua. We see that (p — > ±77 corresponds to a; — s- ±00. 

In the transition region there is some additional energy proportional to the area of the wall. The stress-energy 
tensor is defined as 

T^u ^ id^(j)^)id,<j)^) ~ g^.C (2.10) 



and inserting the solution (2.9 1 we get 



1 , /I , _A , X 



T^u^^^\ri^cosh'\^)di&g{l,{),-l,-l). (2.11) 



Energy and stress of the domain wall are given by 



1 i X 

Pw = {Too) = -A?]'' cosh" ( — ) , 



(2.12) 



That is, the equation of state is 



2 

Pw = -T,Pw (2-13) 



The surface energy density of the domain wall follows as 

+00 

Toodx ^ '^V2X (2.14) 

It is clear that the generation of domain walls with energy density a originates from the spontaneously broken discrete 
symmetry. 

Now we consider a macroscopic volume containing a large number of randomly oriented wall surfaces. For simplicity 
we consider a flat Universe with Robertson- Walker metric ds^ = — R{t)^{dx^ + dy^ + dz^) with scale factor R{t). 
From the Einstein equation we can derive 

with G the gravitational constant. Integrating this equation we arrive at 

R{t)^t\ y^l^Gt^ (2.16) 

This means that the energy density of the walls will quickly overpower the radiation contribution, causing a period 
of power-law inflation. An expansion period of this form would leave less time for galaxy formation. Moreover the 
production rates during nucleosynthesis would change. In case this additional energy density is observable in the 
present Universe, it would cause distortions in the cosmic microwave background, violating the limits on homogeneity 
and isotropy. Quantitative studies show that only surface energy densities of a few MeV are allowed [53] . 



If we trust this cosmological argument we see that the discrete symmetry we encounter in the NMSSM must 
be broken explicitly in order to avoid the formation of domain walls. One idea is to introduce additional higher- 
order operators, that is non-renormalizable breaking terms in the superpotential. These terms are adjusted to 
be Planck-mass suppressed such that they do not have any effect at the low energy scale and renormalizability is of 
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no concern. Nevertheless it can be shown that these additional terms generate quadratically divergent tadpoles for 
the singlet [57t [58] . The generic contribution of higher-order operators to the effective potential with cutoff at the 
Planck-mass scale Mp reads [5T] 



with soft supersymmetry breaking mass parameter rus and ^ depending on the loop order at which the tadpole is 
generated. This contribution to the potential changes the vacuum-expectation-value of the singlet to values far above 
the electroweak scale, that is, destabilizes the gauge hierarchy. 

However one may impose new discrete symmetries in a way to forbid or at least loop suppress the dangerous tadpole 
contributions which arise from the additional Planck-suppressed terms in the superpotential [STJ EHl EHl EI] • Imposing 
a ^2 R-symmetry on the non-renormalizable operators under which all superfields flip sign, the dangerous operators 
are forbidden and there remains only a harmless tadpole contribution to the potential, 



This term breaks the Z3 symmetry and thus avoids the formation of domain walls. In the following we assume that 
the domain wall problem is circumvented in this way without any modifications except far beyond the electroweak 
scale. 



5V = ^MpmlS + c.c. 



(2.17) 



5V = ^mlS + c.c. 



(2.18) 
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3. THE HIGGS-BOSON SECTOR OF THE NMSSM 



In this section we investigate the Higgs-boson sector of the NMSSM, starting with the Higgs potential, electroweak 
symmetry breaking, followed by the derivation of the Higgs-mass matrices and eventually discuss studies on Higgs- 
boson phenomenology. 



3.1. The Higgs potential 



The superpotential (2.1 1 is a holomorphic function which determines all non-gauge interactions of the chiral super- 
multiplets as is discussed in more detail in App. [Cj Here we want to focus on the Higgs sector and first derive the 
physical Higgs potential V and in a second step derive the Higgs mass matrices. With view on Tab. |l] we have the 
scalar fields in the Higgs-boson sector 



Hi] 



s 



(3.1) 



The Higgs-boson potential gets contributions from the chiral supermultiplets, the so-called F-terms encoded in the 
superpotential, from so-called D-terms of the gauge multiplets as well as from the soft supersymmetry breaking terms 
(see App. ICl. 



We start with the scalar part of the superpotential, which reads 



(3.2) 



where we just replaced each supermultiplet (denoted by a hat) in (2.11 by its scalar component according to Tab. |T] 
Here the term (H^eHd) written with the isoweak indices a,f3 — 1,2 reads {Hu)a(-"'^{Hd)[j — H^Hj^ — H^H^, where 
we defined 



1 

-1 



(3.3) 



The so-called F-terms of the Higgs potential are derived from the superpotential in the following way: 



5W 



IAH^P (HlHu + HlHd) + \\{HleHd) + nS''\^ , 



(3.4) 



where 

^=(H„,iJd,5) (3.5) 

denotes the tuple of the three Higgs-boson fields. 

The NMSSM Higgs singlet S is supposed only to couple through the superpotential W to the other Higgs-bosons 
and has no gauge couplings. This restricts the D-terms to be exactly the same as in the MSSM: 

= ^ ^5^(0lT>,)(0jT>,) = \gl\HlHa? +\{9l+ qD {hIh^ - hIh^)^ (3.6) 

id 

Here, T" are the generators of the gauge groups and the corresponding gauge couplings, that is for the group 
U{1)y we have goT° = giYw with the gi the corresponding coupling and the hypercharge operator Yw as well as 
for the group SU{2)l we have ^qT" = g2iw — 52f'^/2, a ~ 1,2,3 with g2 the corresponding gauge coupling. Note 
that we used the identity (jfjcr'^i = "^S.uSjk — SijSki- 

Since there is no knowledge of the realization of the supersymmetry-breaking mechanism (for a discussion of soft 
supersymmetry breaking scenarios we refer to the general literature of supersymmetry and the MSSM), simply all 
possible additional terms are imposed which have couplings of mass dimension one or higher. Only terms which 
violate matter parity are discarded; we briefly recall the terminus of matter parity in App. |C] 
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In the Higgs-boson sector the soft-breaking terms corresponding to the superpotential W are 

f 1 \ 

The first hne is the generic expression, where we get for the scalar fields corresponding mass terms, a general fi-term 
generates corresponding 5-terms, and a trilinear term in W generates corresponding trilinear A-parameter terms. 
Note that there is implicit summation over the indices i,j respectively k. The terms proportional to 6*-' in the first 



line of (3.7) are absent in the NMSSM since we have no corresponding /j,-terms in the superpotential (3.2). However, 
we get the trilinear soft supersymmetry breaking ^-parameter terms, A\ and in the NMSSM. The couplings XA\ 
and KAf^ may be chosen to be real and positive since any complex phase may be absorbed by a global redefinition of 
the scalar Higgs fields 0. Since A and k are in general complex this means that also Ax and are complex in general. 

Eventually we arrive at the scalar Higgs potential of the NMSSM 
V =Vp + Vd + l^soft 

= \X\''\S\''(HiHu + HlHd) + \{HleHd) + KS''\^ 

+ IgilHlHdl' + lial + 9l) (hIHu HlHd)" (^'^^ 

+ ml^HlHu + rnl^HlHd + ml\S\^ + (^XAx{H:^ eHd)S + ^i^A^S^ + c, 

3.2. Tadpole conditions 

We can always parameterize the complex entries of the fields Hd, and S in the the form 

^^=(,^ j' ^" = ^ (,;^K + /.„ + *a.)j' ^==^e^=(.. + /^. + ^a.), (3.9) 

where we can choose the vacuum- expectation values Vd, Vu and Vg to be real and non-negative (any complex phase 
of Vd can be rotated away by gauge transformations and obviously the phases 0„ and <f>g can be chosen to make Vu 
and Vg real and non-negative). Note the convention with a factor l/\/2 in front of the vacuum-expectation- values. 
In the literature this is not handled consistently. Of course the values Vd, and Vg will take the values for which 



the potential V in (3.8) has a global minimum, justifying the terminus vacuum-expectation values. Note that in the 



parameterization (3.9) the vacuum denoted by () precisely means 



{Hd) = ['"^/^ , (Hu) = e'^" {^.J^) ' = e^^'^Vg/V2 . (3.10) 

The stationarity condition for the scalar Higgs potential V at the vacuum translates immediately into the following 
system of equations, so-called tadpole conditions; see for instance |62j . 

9V\ 2 n / 2 2^ ,2 2^ ^ 2 n 

— ) = WdTOff, - R\VuVg + ^ Vd{Vd - W„) + ^fd(w„ + t'J - -^VuVg = , 

dV \ 2 9l + .92 / 2 2\ I'^P / 2 2\ R 2 

- VurriH^ ~ RWdVg + ■ — 5— i'«(i'd - vj + — w«(w<j + w J - i^VdVg = , 



dhu/ ^" 8 2 2 

— \ = Vgml - R\VdVu + RkvI + -^Vg{vl + vl) + [k^w^ - RvdV^Vg = , 

d^J-I^-^^s+^VdVg^O 



I 2 

hvuVg + T;VuVg = , 



(3.11) 



= hVdVu - hVg - IVdVuVg = , 
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where the following abbreviations are introduced 

R = Rc( XK*e'^'^-~^'^^A , / = Im f XK*e'^^--^^^^ 



Rx = ^Re ( A^Ae'^^"+^=^') , h = ^Im ( XAxe''^*-+^' 



^/2 V ^ y ^ V2 



(3.12) 



V2 ^ " 

Note that only these combinations of phases occur in the tadpole conditions. Obviously the last three conditions 



of (3.11 ) can be written in the form 



Ix^~-vs, I^ = -^I (3.13) 

2 2 

and only one of the three imaginary parts /, I\, is not fixed by the others through the vacuum conditions. 

The MSSM may be obtained by taking real parameters A, k in the limit A, k — > with the ratio k/A kept constant 
and the product Xvg as well as the parameters A\ and fixed. In particular this means — > oo in this limit. This 
ensures that the singlet S decouples completely and we end up with the MSSM superpotential. 

3.3. Higgs-boson mass matrices in the NMSSM 

The mass squared matrix of the neutral Higgs scalars is obtained by the second derivative of the Higgs potential 
with respect to the fields at the vacuum. The mass squared part of the corresponding Lagrangian is 

CSSS- - -icV-) i^Zi^ t!:) (t) . (^) . (3.14) 

in the basis — [hd, h^, hg) and ~ (a^, a„, as). The square block-matrices are explicitly 

m|c| + {Rx + Rvs/2)vst(3 |ApWrfU„ - m%spCf3 ~ R^Vs - Rv1/2 \X\'^VdVs - Rxv^ - RvuVs 

Ms = ( \X\'^VdVu - m%spci3 - RxVs - Rvl/2 m%sl + [Rx + Rvs/2)vs/tp lApw.^w^ - RxVd - RvdVs 

\X\'^VdVs - RxVu - RVuVs \XfvuVs - RxVd - RVdVs 2\k\^vI + RxVdVu/Vs + Rf^Vs 

(3.15) 



\Rx + Rvs/2)vst[3 {Rx + Rvs/2)vs {Rx~Rvs)vu 
Mp =\ {Rx + Rvs/2)vs {Rx + Rvs/2)vjtp {Rx - Rvs)vd | , (3.16) 

{Rx-Rvs)vu {Rx-Rvs)vd RxVdVu/vs-3Ri^Vs + 2RvdVu, 



-h'uVs^ 



Msp^I-\ -ivdVs\. (3.17) 

hvuVs hvdVs 2vdVn 



Here the mass parameters m^, m^, mj as well as Ix and are replaced employing the tadpole conditions (3.111. 
From the mass squared matrix we see that in general the neutral Higgs bosons mix. In the case 7 = the mass matrix 
becomes block-diagonal and the scalar fields {hd,hu,hs) do not mix with the pseudo-scalar fields (ad,a„,as), that is 
there is no CP violation in this case in the Higgs-boson sector. Only in this case the neutral Higgs bosons can be 
assigned corresponding CP properties. 

We may isolate the massless Goldstone field G by a change of basis with the /3 rotation (c^ = cos(/3), = sin(/3)) 



^a,d\ ( Cf} S/5 



a,, = 



a J V 1 




Sfj cp Q\ \ a \ = R{/3) a'. (3.18) 



With this rotation, the neutral scalar mass part of the Lagrangian (3.141 reads 



^Higgs mass _ ^(h'^r,'T^\( ^'^^ ^^'sp\(^''\ ( ■i T Q\ 
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with changed matrices 



M'p 





R{PfMpR{f3) = ( (2i?A + Rvs)vs/sm{2(3) {Rx ^ 

(Rx - Rv,)v Rxv^/vs ■ sin(2/3)/2 



Rvs)v 

- SRkV. 



Rv^sin(2P) 



(3.20) 



and 



M'sp = MspRiP) = 



Iv 



'0 -SvsSfj 

—3VsCj3 

,0 Vs ivspcp 



(3.21) 



Since the massless 



Here we introduced the convenient and usual abbreviations tp = tan(/3) = Vu/vd and = 4 
Goldstone boson is now separated (giving the longitudinal mode of the Z-boson) we can suppress the corresponding 
vanishing forth row and column in the mass matrix in ( 3.19[ ). We end up with five physical neutral Higgs-boson 
fields. The mass eigenstates of these five physical fields are obtained by another orthogonal rotation R in the basis 
4> = {hd, hs,a, Qs)'^ yielding 



Hi — R 



with dia,g{m%^,m'jj^ 



'■H31' 



il,,mj,)^RM'R^, 



that is by diagonalizing the mass matrix 



M' 



( Ms 

[{M'spf M'p 



M'sp 



(3.22) 



(3.23) 



In the case of no mixing of the scalar with the pseudoscalar Higgs bosons, the matrix M' is block diagonal and has 
non- vanishing entries only in the upper left 3x3 block and the lower right 2x2 block, which may in this case be diag- 
onalized separately. Then we have three scalar mass eigenstates, denoted by Hi, H2, H^ as well as two pseudoscalar 
mass eigenstates, called Ai and ^2- Without loss of generality the Higgs bosons are finally put in ascending order, 
that is mu^ < T^i/fj < ... < ^tt-Hs or in case of a CP-conserving Higgs-boson sector mu^ < TO_f/2 ^ friHi and toai < toa2- 



Note that in case of k = we get from ( 3.12 ) i? = / = as well as i?^ — — 0. In this case the mass matrix ( 3.23 1 



becomes block diagonal, since M'gp in (3.211 vanishes. Moreover, in this CP conserving case the pseudoscalar mass 
matrix (3.201 has an additional vanishing eigenvalue, giving the so-called axion. This makes clear that a vanishing 



parameter n is accompanied by an additional massless pseudoscalar state. 



The mass squared matrix of the charged fields is defined quite analogously from the Higgs potential with respect 
to the complex (charged) fields, that is in the basis ((i/^)*,77+) the mass matrix reads 



M± = {Rvl/2 + Rxvs + spcpml^ - ]:V^Vd\\\^) - ^^'"''^^^ 



1 



\^ 1 cot(/3) 



Again, the (3 rotation allows us to separate the Goldstone mode 

{H-Y\ ^( C0 s,\ (G+ 
H+ J \-Sf} cp) {h+ 



(3.24) 



(3.25) 



with = {G^y* and H^ = (iJ+)*. From the diagonalization R{j3)'^ M±R{j3) we can read the mass squares of the 
charged Higgs bosons, that is we have two charged Goldstone modes as well as the charged Higgs bosons with 
mass 



4± =tr(Af±) = m^- |ApwV2 + (2i?A + i?i'.)i's/sin(2/3). 



(3.26) 



Note, that by means of this equation we may express the phase parameter Rx in terms of the charged Higgs boson 
mass. For instance, the pseudoscalar mass squared matrix (3.201 reads, using (3.261 



^'^^ ~ \ im|,^ sin(2/3) - \ml^ sin\2P) + sin(2/?) - iR^i 



(3.27) 



with the abbreviation for the upper left entry 



(3.28) 
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bosons 


gauge eigenstates 


mass eigenstates 


sleptons 


Cl, Sr, i>e 
jj-L, P-R, ^fi 
TL, Tr, Vr 


Cl, Cr, 
P-L, P-R, t'n 
Ti, r2, Ut 


squarks 


UL, Ur, dL, dR 
Cl, Cr, Sl, Sr 

iL, tR, Bl, bR 


UL, UR, dL, dR 
Cl, Cr, Sl, Sr 

ti, t2, bl, 62 


Higgs bosons 


hd, hu, hs, CL, CLs 

Hd , Hu 


Hi , H2 , Hi , H4 , Hs 
{Hi, H2, Hi, Ai, A2) 
H^ 


fermions 






neutralinos 


B°, W\ Hi, H'l S 


X?, X2, X3, X4, xi 


charginos 


W^, Hd, H+ 


xt, X2 


gluino 


9 


9 



TABLE II: New particles in the NMSSM compared to the SM. The gauge eigenstates as well as the through mixing generated 
mass eigenstates are shown. The Higgs-boson mass eigenstates in brackets give the notation in the case of a CP-conserving 
Higgs-boson sector, where we have three CP-even Higgs bosons Hi, H2, Hi as well as two CP-odd ones, Ai and A2 instead of 
five mixed states Hi, ...,Hs- 



In the MSSM limit of the model, we have no mixing with a singlet Higgs boson state as and the upper left 
mass squared entry becomes the single pseudoscalar mass squared m?p = (m^^^^)^. But in general even in the 
CP-conserving Higgs case, the two pseudoscalar mass eigenstates in the NMSSM, niAi and niA^, originate from the 



mixing in (3.27 1 



The only modifications of the NMSSM compared to the MSSM arise from the extended superpotential, that is the 
generalized /i-term, imposing a gauge singlet. Nevertheless, the mixing of the corresponding additional states with 
other states may generate a coupling of the singlet to the gauge bosons. In table |IT] all new particles predicted by the 
NMSSM are shown. The second column lists the gauge eigenstates whereas the right column gives the corresponding 
mass eigenstates. 



3.4. Stability and electroweak symmetry breaking of the global minimum 



In this section we want to present conditions of the parameters in the NMSSM which follow from stability, as well 
as from the required symmetry breaking behavior at the vacuum. First of all, the Higgs potential has to be bounded 
from below in order to yield a stable vacuum solution. From the F-terms as well as from the D-terms of the potential 
we get quartic terms in the Hd and Hu fields with positive coefficients independent of the choice of the A parameter. 
These quartic terms dominate the potential value for large field values. The singlet S has only one quartic F-term 
with positive coefficient as long as k, is non-vanishing. That is, stability is guaranteed in the NMSSM potential if 
K ^ 0. 

Of course the potential has to have a minimum with the right electroweak symmetry breaking behavior. In particular 
this means that the potential value at the vacuum has to be lower than at the symmetric stationary point with 
{Hd) = {Hu) = 0. Moreover a stationary point with {S) = has to be avoided in order to generate an effective /i 
term and thus gives masses to up- and down-type fermions. Since the Higgs potential is zero for vanishing scalar 



doublet and singlet fields we must have {y) < 0. Inserting the parameterization (3.9) into the potential (3.8 1, using 



the tadpole conditions and the replacement of R\ by the charged Higgs mass (3.26) we get an upper bound 



m 



2 r, 1^1 9 l^l'^'^t 4i?KWs 2 

^± sinM + «2 sin(2/3) + tan2(2/3) ~ sinM ^ sin^m ^ """^ ' 



Note that this upper bound goes to infinity in the limit — s- cx), that is in the MSSM limit there is no upper bound. 

A further constraint of a minimum of the potential is to have a positive definite Hessian matrix in the scalar fields, 
that is a positive definite scalar mass-squared matrix at the vacuum. Since the eigenvalues of the Hessian matrix are the 
mass eigenstates, they have to be positive. In the simplified CP-conserving case where the scalar mass-squared matrix 
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becomes block diagonal, this condition requires det Mp > 0, that is (rn^± — m^+l/2|ApWs)(3/4i?t;^ sin(2/3) — > 
?>/4:E?v'^vl, giving 



(3.30) 



That is in the NMSSM we find for non-vanishing A a reduced lower bound for the charged Higgs-boson masses 
compared to the MSSM. The conditions (3.291, (3.301 together with the tadpole equations are necessary conditions 



for a viable global minimum, but not sufficient. A stringent determination of the global minimum which has the 
required electroweak symmetry breaking behavior can be found in (53] and will be discussed in more detail in Sect. [6] 



3.5. Parameters of the NMSSM Higgs potential 



The parameters of the NMSSM Higgs potential (3.8 1 are 



(3.31) 



in addition to the SU (2)^ x U{1)y couplings 52 and gi . From reality of the Lagrangian we see that only A, k, Ax, A„ 
are complex whereas Tnjj^^rrijj I'^^s have to be real. In practice, it is often not useful to fix these initial parameters, 
but others like the vacuum expectation values of the Higgs fields, Vd, 

parameters , m'^ 

the diagonalization of the mixing matrices as discussed in Sect. 3.3 



or typically, written via tan (3 



for the two Higgs doublets in terms of v 



246 GeV, tp, Vg- Moreover, we remark that the 
are not the physical parameters of the Higgs bosons since the physical masses arise from 



The tadpole conditions (3.11) and relation (3.26) can be employed to translate the new set of parameters, 



A, K, \Ai^\, V, tan (3, Vs, m^±,signi?K, S] 



EDM, S'^ 



(3.32) 



into the initial parameter set (3.31). Here the phase combinations 



JEDM 



(3.33) 



We 



are defined and the complex parameters are writte n as A — |A|e^^, k — \K\e " , A\ = \A\\e ^^^j A^ = \A^\e 
can easily see how to get the original parameters (3.31) back from this new set: via the phase combinations ^edm 
and (5^ we can determine R and / in (3.12). The tadpole conditions fix then I\ and 1^. Together with the l ength 
the phase parameter i?^; is determined. The physical charged Higgs mass 'rrc\j± gives via (3.26) 

Together with u, tan /?, Vs we get all 
lie conditions. It only remains to fix the phases of A\ and A 



l^^l and signi?K 

immediately that is all phase parameters (|3.12|) are known at this point, 
mass parameters "tv?^^ , rri}^^ , m| from the tadpo 



\AJ\ are known. Note, 



to arrive 
that in the mass matrix 



at the initial set. This is done by using again (3.26 1 since |A|, |k 
of the Higgs scalars, the CP- violating entries in the matrix M'^p in (3.21) are proportional to the imaginary part of 
exp[i(5EDM - (^k)]- 



3.6. The one-loop effective potential 



For later convenience let us also introduce the contribution of radiative corrections of the Higgs-boson mass to the 
potential, which may be in-cooperated by studying the effective potential 

Kff - + Moop . (3.34) 

The contributions of corrections from the light particles, that is leptons and the quarks and squarks of the first two 
families may be neglected due to the small Yukawa couplings to the Higgs fields. At the one-loop order, the result for 
Vioop, first given by Coleman and Weinberg [53], reads 

^oop = ^ E (-1)'^* (25. + 1) < In (^) , (3.35) 

i 

where the sum is over all particles with field-dependent masses m^, spin Si and color degrees of freedom Ci and Q is 
the renormalization scale from the loop integral. Of course, this corrections of the Higgs potential change the tadpole 
conditions. The field-dependent masses, that is the masses before spontaneous symmetry breaking occurs, of the top, 
the stop and the gauge bosons are given in App. |F] 
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3.7. Higgs-boson phenomenology 

In contrast to the minimal supersymmetric extension of the SM we have in the next-to- minimal extension, NMSSM, 
due to the additional singlet superfield S two additional Higgs bosons as well as one additional singlino. In particular 
the phenomenology of the Higgs bosons may be very different compared to the MSSM. This difference is pointed out 
in this subsection. We start with considering one of the advantages of the NMSSM over the MSSM, namely the much 
less restricted Higgs-boson sector. In the MSSM the lightest scalar Higgs boson is predicted to not exceed the Z-boson 
mass, 

(m^f™)2 <m|cos2(2/?). (3.36) 

This bound is softened through quantum corrections. One-loop corrections to the lightest Higgs boson where stud- 
ied [SH |66l Ell [68] as well as dominant two-loop corrections [39l [68] • The largest contribution comes from virtual top 
and stop loops, where the one-loop corrections read 

(Am^f SM)2 ^ In f !!%!^) , (3.37) 

with c of the order 1. Confronted with the experimental data, that is LEP searches for light Higgs bosons in the 
MSSM we have the lower bound [69] 

^MSSM. exp ^ 92 GeV . (3.38) 



That is, we need, considering (3.361 and (3.371 rather large stop masses, which enter only logarithmically in 
the radiative corrections, in order to increase the mass of the lightest scalar Higgs-boson in the MSSM. On the 
other hand, a large violation of the degeneracy between the superpartners top- and stop-mass reintroduces a 
new kind of unnatural fine-tuning. This is in particular disturbing since one of the main motivation for the 
introduction of supersymmetry is to avoid this unnatural large quantum corrections to the scalar Higgs-boson 
mass, which occur in the SM. Let us note that there is some debate about this argument and refer the reader to 
the discussion in Martin [2T]. We notice, that a large mass splitting of the superpartners would reintroduce the 
initial naturalness problem we encounter in the SM. Therfore, we expect the superpartner particle masses not to 
exceed the scale of 1 TeV too much to be considered as natural. Another argument for superpartner masses not 
to exceed the scale of 1 TeV to much is the unification of gauge couplings, which is spoiled by very large mass splittings. 



The situation in the NMSSM is quite different. Firstly, the tree-level mass bound (3.36 1 for the lightest Higgs boson 



is no longer valid in the NMSSM. In the CP-conserving Higgs-boson sector case the lower bound is changed to [70] 

(^NMSSM)2 < ^1 / 2(2^3) + ^-^f^) . (3.39) 

\ yi + .92 / 



Thus, the upper limit of the lightest CP-even Higgs-boson mass is in general lifted compared to the MSSM (3.361 



Secondly, the experimental bound ( 3.38 ) is also much weaker in the NMSSM [7TJ [72J [73] . This can be easily understood: 
The main detection strategy of the lightest scalar Higgs in the MSSM is based on Higgs-strahlung off a Z-boson with 
subsequent decay of the Higgs boson according to Hi bb and fe-tagging in the detector; see Fig. [l] But in the 




FIG. 1: Higgs-boson production via Higgs-strahlung and subsequent decay into a bb quark pair, as searched for at the LEP 
experiment '691. The coupling of the Htggs boson to the Z bosons is denoted by qhzz and given for the NMSSM m App. [S[ 

NMSSM the additional decay channel into a pair of pseudoscalar Higgs bosons. Hi — + AiAi, is open, reducing in 
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general the bb decay-channel branching fraction. In the case m^i < 2mf,, the subsequent decays of the pseudoscalar 
Higgs bosons can no longer proceed via b-quark pairs. In this case, the pseudoscalar decays into t^t~ , cc or gg and 
may escape detection. Then only the weaker decay mode independent LEP bound [71] 

^(N)MSSM, cxp ^ 32 GeV (3.40) 

applies. This experimental search is not based on the decay products of the Higgs boson, but on the recoil mass 
spectrum of the Z boson in the process e+e^ ZHi. Therefore, the actual decay of the Higgs boson plays no 
role in this search and the given limit applies also to the NMSSM. To summarize, in the NMSSM we have a larger 
theoretical upper bound as well as a lower experimental bound and thus this model seems to be favored over the 
MSSM with respect to the Higgs-boson sector restrictions. 

With these remarks on the restrictions in the Higgs-boson sector, let us also briefly comment on the approach of 
Barate et al., [75], on Higgs-boson phenomenology in extensions of the SM. In Fig. [2] the upper limit on the ratio 

times the branching ratio BR(i/ — > bb) is shown, depending on the Higgs-boson mass. In the ratio the expression 
5hzz denotes the SM Higgs-Z-Z coupling, whereas gnzz denotes the same non-standard coupling. The dark and 
bright bands give the l-cr and 2-(t deviations from this limit. The limits on are gained from LEPl data at the 
Z-resonance as well as from LEP2 data taken at energies between 161 and 209 GeV. The Higgs boson is assumed 
to decay into fermions and bosons according to the SM. The dominant decay channel in the SM proceeds via the 
bb channel. The Higgs-boson production cross sections for the processes e^e" HZ, WW H and ZZ H 
are scaled with the non-standard coupling squared g^2,z- Obviously the observation exceeds the upper limit for an 
assumed Higgs-boson mass of around 98 GeV. As is argued in the works of Dermisek et al. [7TJ |72], the observed 




FIG. 2: Upper limit on the ratio ^ given in ( |3.41 1 times the branching ratio of the Higgs boson into bb according to the SM. 
The upper limit is gained from at LEP from LEPl data at the Z resonance as well as from LEP2 data taken at energies between 
161 and 209 GeV. The dark and bright band give the l-o and 2-a deviations from the central upper limit. The observation 
is shown, exceeding the upper limit for a Higgs-boson mass of about 98 GeV. The horizontal line above represents the SM 
expectation on . Figure taken from \7S^. 



excess is consistent with models which have a SM-like ZZH coupling but a reduced branching ratio BR(77i — > bb) 
due to new open decay channels. As aforementioned, in the NMSSM the light scalar Higgs boson may decay into 
a pair of pseudoscalars Hi — > AiAi. In case niA-^ < 2mi, the pseudoscalar Ai can not subsequently decay into bb 
pairs, but decays into r pairs or jets. As is argued, for a SM-like ZZH coupling and BK{Hi bb) « 0.08 and 
BR(i/i AiAi) « 0.9 the observed excess at rriHi « 98 GeV is nicely reproduced. 



Let us turn our attention to the prospects of NMSSM Higgs-bosons at the LHC, where we mention the approaches 

[33 [Ml [Ml IZi [7i IZZl IZHl El [Sni ini [HI [HI ^ ^ ,8^ • 
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In Ref. [75], published in 1996, the Higgs detection capabiHty of LEP2 and LHC was studied. It was stated, that 
there is large parameter space, where all Higgs bosons in the NMSSM escape detection at both collider experiments, 
LEP2 and LHC. With respect to the LHC, an integrated luminosity as high as 600 fb~^ was assumed. 



However, in the meantime much improvements compared to the initial study could be achieved. In 2002 an 
investigation of Ellwanger, Gunion and Hugonie [83], was titled rather promising: Establishing a No-Lose Theorem 
for NMSSM Higgs Boson Discovery at the LHC. At this time there were improvements at the theoretical side, like 
the two-loop corrections to the effective potential [321 IHH] and the predictions of the VF-fusion channel Higgs-boson 
production at the LHC [TT] |7H1 [ZSl EHl EI]- Moreover the improved LEP2 data were available, constraining the 
NMSSM parameter space further. Taking the newly predicted detection channels into account, that is the associated 
production of Higgs bosons with a top pair and the two VF-fusion channel, altogether the following channel list is 
studied in |83j : 

1) 5.9 ^ -ff ^ 77 , 

2) WH or ttH production with 77Z^in the final state , 

3) tin with jH ^bb, 

4) gg ^ H/A or bbH/A production with H/A tt , 

5) gg^ H ^ ZZ* 4 leptons , (3.42) 

6) gg ^ H ^WW* ^l+rvv , 
at LEP2: e+e" ZH and e+e" HA, 
WW ^ H ^Tf, 
WW H ^ WW* , 

where H denotes any CP-even and A any CP-odd Higgs-boson, that is a CP-conserving Higgs-boson sector is 
considered. In this study, the Higgs-to-Higgs as well as Higgs-to-top decays are not taken into account, that 
is parameter sets leading to a particle spectrum which allows kinematically for these decays are disregarded. 
Exphcitely, the disregarded decay channels are H -> HH, H AA, H H+H~, H AZ, A HA, A — > HZ, 
H/A H^W^, H/A ^ tt, t ^ H^b. The suppression of these decay channels is implemented by invoking the 
following constraints: 



niHs < 2mHi,2mAi,2mH±,mA:^ + mz,mH± + mw, 

< niHi + rriAi , rriHi + mz, mH± + mw, 
mH± > 155 GeV. 



(3.43) 



Further assumptions are the absence of Landau singularities (see Sect. 5.11 and, as mentioned before, the LEP2 



constraints on Higgs strahlung e'^e ZH [87 



HA 



Also in the parameter space it is required that 



1^1 = Xvs > 100 GeV in order to avoid a light chargino, which is experimentally excluded. Neglecting kinematically 
the Higgs-to-Higgs, Higgs-to-top decays and Higgs-to-neutralino decays, it is found that at least one Higgs boson 
will be detected for a luminosity of 300 fb^^ for arbitrary choices of remaining parameters, with a statistical discovery 
level of at least 5-(t. 



In 2003 this investigation was followed by supplementing the missing H AA decay channel via VF-fusion, that 
is WW —^H—f AA |84l (Towards a no-lose theorem for NMSSM Higgs discovery at the LHC). There, special 
parameter sets are chosen, representing cases which are rather difficult for detection at the LHC in this channel due 
to difficulties with respect to the background. 

In 2005 follows a study |89j {Difficult scenarios for NMSSM Higgs discovery at the LHC) with an investigation 
aiming to reveal parameter sets which yield a mass and coupling spectrum such that no Higgs boson is detectable at 
the LHC. A CP-conserving Higgs sector is considered. Gaugino mass unification is assumed, where M2 = 1 TeV is 
fixed at the electroweak scale, corresponding to Mi w 500 GeV and Mg « 3 TeV. For the soft supersymmetry breaking 
parameters it is set mq = rriu = = rriL = rrie — 1 TeV for all three generations. Further, the A-parameters are 
fixed to Au = Aj^ = 1.5 TeV for all generations. Eventually, a minimal charged Higgs-boson mass, rnH± > 155 GeV, 
is assumed in order to avoid detection of charged Higgs bosons in decays like t H'^b for moderate values of tan(/3). 
This choice of parameters suggests that the Higgs-boson detection might be the only new signal at the LHC. Over 
the following ranges of remaining parameters is scanned: 

10-4 <A< 0.75, -0.65 < At < 0.65, 1.6 < tan(/3) < 54 , -1 TeV < /x = Aw^, ^a, ^« < 1 TeV . (3.44) 
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P 



P 




FIG. 3: Central exclusive production of a light CP-even Higgs boson (H) in proton (p) collisions, pp — > ppH, as suggested 
by 180/ . The Higgs boson H decays into a pair of pseudo scalar Higgs bosons A. The process proceeds via a fermion loop. 



In this study the "no-lose" theorem is confirmed, if the Higgs-to-Higgs decays are kinematically excluded via (3.431. 
As is noted, this exclusion corresponds to large parts of available parameter space. In the complementary region of 
parameter space in deed parameters are found, where no detection at LHC may be achievable, as is reported. We 
repeat two example parameter sets, denoted by set 7 and 8, 



parameter set 


A 


K 


tan(/3) 








7 


0.5 


-0.15 


3.5 


200 


780 


230 


8 


0.27 


0.15 


2.9 


-753 


312 


8.4 



which might lead to a particle and coupling spectrum not detectable at the LHC. In the parameter set 7 the 
next-to- lightest CP-even Higgs boson H2 has SM-like couplings and decays dominantly via H2 HiHi, with 
branching ratio BK{H2 HiHi) = 0.93. Moreover, the Hi has only small couplings to bb and r+T~, following from 
the fact that Hi is mostly iJ°-like in the mixing matrix. This means that the signature WW H2 ^ 2j t^t^ 
is suppressed and there remains only the 4j final state which is of course heavily plagued by QCD background. In 
parameter set 8 the lightest CP-even Higgs boson Hi is SM-like and decays exclusively via Hi AiAi. In this 
scenario the pseudoscalar Ai Higgs boson is very light, that is niAi = 1 GeV, and decays only into light jets. Thus, 
also this signature is overpowered by QCD background. The authors note that at a complementary e+e^ collider, 
like the proposed ILC, the scenarios, difficult to detect at the LHC, could be observed in the decay-independent 
signature e'^e~ ^ — > ZX, as discussed in the end of this section. 

In the work |86| (Reinstating the 'no-lose ' theorem for NMSSM Higgs discovery at the LHC) a clear statement of a 
5-(T discovery of at least one Higgs boson is given. This statement is achieved by employing an additional constraint: it 
is argued that absence of large fine-tuning in the NMSSM corresponds to parameter space with a CP-even Higgs boson 
in the mass range 90 GeV < ttlh < 100 GeV. In this work fine-tuning is studied on a quantitative basis, introduced 
later in Sect. |5.1[ It is argued, that for a CP-even Higgs-boson with mass around 98 GeV, the only remaining decay 
not excluded by LEP may proceed via a light pair of pseudoscalars with subsequent decay into r's or jets, 

H ^ AA^ T+T^T+T- or 4j . (3.45) 

The current experimental lower limits of this process with the decay into final r's is jtj^'^ssm, oxp ^ gg GeV [55] 
and with the decay into final state jets is jt^^^^^m, oxp ^ g2 QeV [71] from the decay mode independent searches. 



This means that Higgs production with subsequent decay (3.451 is well beyond the constraints at LEP. Further, it is 



noted that the requirement < 2mf, is natural in the sense of low fine-tuning. Nevertheless, it is found that this 
decay channel is very challenging to detect at the LHC in all so far studied Higgs production processes. However 
as an additional search strategy, the authors suggest to consider the well-known central exclusive production of a 
light Higgs boson H via the process pp ppH. This process with the subsequent decay into a pair of pseudoscalar 
Higgs-bosons is shown in Fig. [3j This detection mode requires the outgoing protons p to be detected in purpose-built 
low-angle detectors (FP420) [90] IM]. Via these special detectors the four- momentum of the central system could be 
reconstructed very accurately and thus the masses of both the H and the A be determined on an event-by-event 
basis. The investigation shows that this technique in deed would allow for a discovery of the lightest Higgs bosons 
even when background is taken into account. 



The decay chain (3.451 with rriAi < ^m^ and large branching fractions BR{Ai r+r ) was also studied in the 



publication of Belyaev et al. [92| in order to close the gap for a 'no-lose' theorem. The parameters of the Higgs-boson 



sector (see Sect. 3.5 1 were varied in large ranges and the remaining parameters, which enter the Higgs-boson 



sector only at the suppressed loop level, were fixed in this study. The authors were looking for signatures, where 
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two r's decay into /x's and the other two t's into jets. The Higgs-strahlung as well as the vector-boson-fusion 
production channels were considered and, as is reported, thousands of events are predicted, mostly originating from 
vector-boson-fusion processes. An extended investigation, taking also the background into account, is announced to 
be in preparation. 

Eventually we note that at a future electron-positron collider like the international linear collider (ILC) with 
center-of-mass energies up to 1 TeV, the study of the Higgs-strahlung process e+e^ ^ Z* ^ ZH would be accessible 
for a large mass range of the Higgs bosons. In the ILC reference design report, referring to an investigation in 
the constrained NMSSM [53] , it is stated that the measurement of the Higgs-boson masses with a resolution of the 
order of 100 MeV could be achieved, if the Higgs bosons are not too heavy. A very direct approach available at an 
electron-positron collider is the Higgs-boson detection, independent on the subsequent decay-mode. From the recoil 
mass spectrum against the Z-boson in the process e+e~ ^ Z* ^ ZX the Higgs boson is detectable if the couplings 
to the Z-boson are not too small. This decay-mode independent search is in particular interesting in the case, where 
the Higgs-to-Higgs decay is dominant and the LHC might fail to detect the decay products. 
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4. NEUTRALINOS 



The modification in the NMSSM compared to the minimal supersymmetric extension originates from the superpo- 
tential terms XS{H^€Hd) + ^^S^ accompanied by the introduction of an additional singlet S. Here we focus on the 
modification which arises from the fermionic part of the singlet S, the so-called singlino S. For studies of the MSSM 
phenomenology we refer the reader to the review of Martin [21] including detailed references. Note that we do not 
have additional charginos compared to the MSSM. The singlino component of the superfield gives a fifth neutralino 
which in general mixes with the bino B, wino W'^ and Higgsinos and see Tab. |ll] Collecting all quadratic 
terms as performed in App. |D]we arrive in the basis ip'^ = {B,W^ , Hj^, H^, S)"^ at the symmetric neutralino mass 
matrix 



Ml 


-cpswniz 
s/jswrnz 






Ah 
cpcwmz 
-spcwmz 




-cpswmz 
ciscwmz 




spswmz 
-spcwmz 
-Xvs/V2 


-Xvd/V2 



' 


y/2KVs 



(4.1) 



This neutralino mass matrix refiects the fact that the singlino S does not couple to the gauge bosons but only 
to the Higgsino doublets H^, in addition to the selfcoupling in the lower diagonal entry. In case of a small 
mixing for the singlino with the Higgsinos the singlino decouples. In this case the behavior of the four neutralinos 
is MSSM- like and the mass of the singlino is at tree level given by the lower diagonal entry m| « 2k^v'^. For large 
mass values mj the singlino may escape detection and a distinction of the NMSSM from the MSSM is very difficult [94] . 



After diagonalization with the unitary rotation U of the mass matrix (4.1 1 (see (D27l for more details) we end up 
with five neutralinos 



with i, j = 1, 5 



(4.2) 

which afterwards are arranged in ascending order, thus we end up in Dirac notation with « = 1, 5, with Xi the 
lightest neutralino. The singlino components of the mixed states do not couple to gauge bosons, gauginos, leptons, 
sleptons, quarks and squarks. Thus, in addition to an increased total number of neutralinos, we expect in particular a 
changed behavior compared to the MSSM for a neutralino with a large singlino component, that is for a singlino-like 
neutralino. It is worthwhile to note, that the detection of a fifth neutralino would be a clear signal of an extended 
supersymmetric model. However, the production cross section of a singlino-like neutralino is rather small due to its 
small couplings. Moreover it is evident from its small couplings that a singlino-like neutralino, which is not the LSP, 
would be omitted in cascade decays. 



In Fig. |4] the mass spectrum of the five neutralinos as well as the dominant mixings (couplings) of the mass 
eigenstates are shown depending on fix = Xv/V2, as presented in Ref. [94 . The parameters chosen for this plot are 
given in the figure caption. The parameter ii\ determines the mixing of the singlino with the Higgsinos, as is evident 
from the mixing matrix (4.1 1. We see that in this scenario the lightest neutralino is singlino-like, since the mixing U^^ 
is dominant for small values of fi\ as shown on the right hand side of Fig. [4] This is what is expected since for small 
fix the singlino decouples from the other neutralinos. Let us note that in Ref. |94j a study of neutralino production 
cross sections in e+e~ collisions and decay rates of the neutralinos can also be found. 



In the study by Hesselbach and Franke [IS] the associated singlino-like neutralino production is discussed 



e^e 



XiXs ■ 



(4.3) 



where denotes a singlino-like neutralino. From the neutralino mixing matrix (4.1 1 we see from the lower diagonal 
entry that we get a neutralino with a large singlino component for a large vacuum-expectation- value Vs ■ As shown in 
Fig. |5] the total cross section drops with increasing singlet vacuum-expectation- value Vg (denoted by x in this work), 
as is evident from the increasing singlino-component in that is smaller couplings for rising Vg- The effect of beam 
polarization is also shown in this figure and we see that we get total cross sections of the order of 10 femtobarn for 
not too large Vg, depending also on the beam polarizations. The chosen parameters in this example are given in the 
figure caption. 



In case four light neutralinos are detected, the distinction of the minimal supersymmetric model from further 
extension is in general difficult. As noted already, a heavy singlino-like neutralino would be omitted in cascade decays 
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FIG. 4: Neutralino masses (left) and mixings (right) devendma on a \ = Au/\/2- The other parameters are fixed to ^finva = 
120 GeV (this is the lower diagonal entry in the mixing matrix 4-lK M\ = 250 GeV, M2 = 500 GeV, \Vs/y/2 = 170 GeV, 
tan(/3) — 3. The numbering i = 1,...,5 denotes the mass eigenstates in the mixing x? = Ujiip^ (4.21, before the arrangement 
in ascending order. The mixing entries m the matrix Uji are denoted in the right plot by the indices, that is for instance |53| 
denotes {Us^l- The dashed curves show an approximation, in case the singlino-Higgsino mixing is suppressed, that is for v ^ Vs- 
Figure taken from \94'l - 




FIG. 5: Total cross section for singlino-like neutralino production e'^e~ XiXs depending on the vacuum- expectation-value 
Vs, here denoted by x as presented by Hesselbach and Franke 195). The electron-positron center- of -mass energy is fixed to 
500 GeV. The model parameters chosen are inspired by the SPSla scenario in the MSSM with Mi — 99 GeV, M2 = 193 GeV, 
tan(/3) — 10, jj, — Xva = 352 GeV. The full line refers to unpolarized beams, the dashed line to beam polarizations P- = +0.8, 
P+ — —0.6 and the dotted line to beam polarizations P_ — —0.8, P+ — +0.6. The mass of the neutralino x's fixed at 120 Ge V 
by the K-parameter. 



of even heavier supersymmetric partner particles. In case there is substantial mixing in the neutralino sector there 
is a method discussed by Choi et al. [31] to discriminate the NMSSM from the MSSM at a future electron-positron 
collider: the total cross section, summed over the four light neutralinos and normalized to its asymptotic form in the 
limit of infinite center-of-mass energies, is sensitive to the total number of neutralinos. The energy dependence of the 
ratios for the MSSM and the NMSSM is shown in Fig. [6] The model parameters chosen in this plot are given in the 
figure caption. 

There remains the possibility that the singlino-like neutralino is the LSP as discussed in more detail in Sect. |5.4[ 
Parameter scans in the NMSSM show that there is indeed large parameter space with this possibility passing all 
nowadays known theoretical and experimental constraints. Due to conserved matter parity this LSP singlino-like 
neutralino would escape detection. However, due to its small couplings the next-to-lightest supersymmetric particle 
(NLSP) would decay very slowly into the singlino-like LSP. This opens the possibility to observe displaced vertices 
[nZl [Ml [nS]- The partial decay width of a sfermion / decaying into the LSP neutralino and a fermion / is jlOOl 
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FIG. 6: Energy dependence of the four light neutrahno production cross sections summed up and normalized to its asymptotic 
form for infinite center- of-mass energies. The full Itne corresponds to the MSSM, whereas the dashed line corresponds to the 
NMSSM. The assumed parameters are tan(/3) = 10, |Mi| = 100.5 GeV, M2 = 190.8 GeV, = \\vs\ = 365.1 GeV where CP- 
invariance in the Higgs-boson sector is assumed. The selectron masses are fixed to rue^ = 208.7 GeV and rng,^ — 144.1 GeV. 
Figure taken from \90f . 

mmnm\ 

r(/ - X?/) - ^^^J; ((a' + b'){mf -m^o-mf)-Aab m^m^o) (4.4) 

with kinematic function y, z) = + + ~ 2xy — 2xz — 2yz and a and b the left- and right couphngs in the 
Lagrangian 

^//x? = + bPL)x"J + c.c. (4.5) 

The length of flight of the sfermion is simply Ij = hc/T{f Xif)- Scenarios of the NMSSM are studied for which 
a NLSP X2 decays into a singlino-like LSP Xi- As is pointed out, rather large vacuum-expectation- values Vs are 
required in order to get a very pure singlino-like neutralino and suppressed couplings accompanied by observable 
displaced vertices. 

In the approach of Kraml, Raklev and White [103' a special scenario is studied, motivated by the SPSla in the 
MSSM and extended to the NMSSM. In this scenario, the LSP is a singlino-like neutralino Xi and the NLSP a bino- 
like neutralino X2 with a small mass difference Am = rri^o — rri^o . Their analysis is based on cascade decays which 
eventually end up with the decay 

X^^X?Z+r. (4.6) 

Due to the small mass difference the final state leptons are soft and could escape detection. As the authors emphasize, 
with a typical kinematical cut on the minimal transversal momentum of leptons at the LHC this could lead to a wrong 
interpretation of a X2 LSP, since non of the true final state particles would be seen in this case. Accepting low momenta 
of the final state leptons, the measurement of the invariant di-lepton mass squared distribution mfi = {p^ -\- p^)"^, 
with the momenta of the leptons denoted with , would allow to determine the mass difference Am = m^o — m^o 
at the edge of a measured wi^i distribution. 
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5. PARAMETER CONSTRAINTS 



In this section we will discuss some parameter constraint studies, performed in the NMSSM. We start with theoretical 
constraints. Constraints originating from a viable global minimum of the Higgs potential, that is a global minimum 
which has the correct electroweak symmetry breaking behavior, are discussed in detail in Sect.|6] We mention the 
requirement of perturbativity of couplings in considering renormalization group equations. Then we turn to fine- 
tuning. In a series of recent papers some attention is payed to fine-tuning in extensions of the SM; for instance in 
Refs. |71l 1721 [75] . where fine-tuning is quantified in terms of a simple function. We proceed with constraints coming 
from the requirement of a vacuum which is color- and electric charge breaking invariant. Also an approximative 
constraint can be gained from the condition to have a non- vanishing vacuum-expectation- value for the singlet (S) — Vg- 
On the experimental side we first consider the firm exclusion limits from collider experiments. Further constraints with 
respect to the anomalous magnetic moment of the muon, the b-meson decay, cosmological constraints from indirect 
as well as direct detection of cold dark matter and strong first order electroweak phase transitions are considered. 
Finally, we review some generic parameter scans over large parameter regions, revealing the viable parameter space. 



5.1. Theoretical constraints 

Let us start with the Landau pole exclusion constraint. The renormalization group equations for the NMSSM have 
been determined to two- loop order. Here we present the one-loop results of the gauge couplings and the superpotential 
parameters A and n as well as yt in the case of a CP-conserving Higgs sector |104j . 

di^'^^ie^y^- = 1,2,3), 



d 2 _ 1 



(5.1) 



(A^ + K^) 



The coefficient for the couplings are hi — 33/5, 62 = 1, ^3 = —3. The context of the gauge couplings gi, (72 with 
the conventional couplings is gi = yjbjigi, §2 = 52 with, as usual, e = g2Sm{9w) = giCos{9w)- The scale factor is 
defined a,s t = ln{Q^ / ttigut) with unification scale togut ~ 3 • 10^^ GeV and Q the scale under consideration. Note 
that the unification of the gauge couplings at the GUT scale occurs also in the NMSSM, like in the MSSM. This is 
due to the fact, that the additional superfield 5 is a gauge-singlet. 

The requirement of perturbativity of the couplings means that quantum corrections are assumed not to become too 
large at all scales from the electroweak scale of about 100 GeV up to the gand unification scale of about 3 • 10^^ GeV 
(typically, the couplings are required not to exceed 2tt). At least there should not be any Landau pole for the 



coupHngs. Due to the /3 functions given in (5.11 we see that A and k drop with a decreasing energy scale t. That 
is, perturbativity up to the GUT scale gives a strong constraint of these coupling parameters at the electroweak 
scale. As an example, in Fig. [7] the running of the couplings A and k with the assumption of real parameter values 
is shown. It is worthwhile to note, that the constraint of perturbativity up to the GUT or Planck scale is based on 
the assumption, that the NMSSM is valid up to the very high GUT or Planck scale. In spite of the unification of 
gauge couplings, which indicate that this could be true, this is not granted even if the NMSSM may be correct as an 
effective theory at the electroweak scale. 

One of the main motivations for the introduction of supersymmetry is to avoid the fine-tuning problem we encounter 
in the SM. That is we expect the model under consideration not to reintroduce fine-tuning again. In order to quantify 
fine-tuning Barbieri, Riccardo and Giudice introduced the quantity |105j . 



F = max 



dlii{m, z{ai)) 



dln(ai) 



(5.2) 



with Qi denoting all soft supersymmetry breaking parameters. Note, that in the original work the Z-boson mass 
occurs squared in this expression. With help of the quantity F, parameter sets within a certain model may be 
compared on a quantitative basis and the parameters giving lower values of F correspond to a more natural parameter 
choice. Moreover this allows also for a comparison of different models with respect to fine-tuning. This approach is 
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Q [GeV] 



FIG. 7: Evolution of the dimensionless NMSSM parameters A and k taken to be real 1104^ - The Yukawa coupling ht is set to 
0.8 at the grand unification scale, Qgut ~ 3 • 10^® GeV. For the evolution of X we start with \{Qgut) ~ 3 and k.{Qgut) ~ 1- 
The evolution of k, is given for vice versa values, that is k{Qgut) = 3 and \{Qgut) = 1. See Ref. fW^ for more details. 



employed in some studies on the Higgs-boson spectrum in the NMSSM; see Sect. 5.4 



A further restriction for the parameters in the potential comes from the requirement of a color and electric charge- 
invariant vacuum [HUl 1106] . This is evident if we consider for instance the scalar part of the slepton-slepton-Higgs su- 
perpotential term —e*pye{L'^tHd). With view on Tab.|l]we see that the supermultiplets have hypercharges Y{eii) — 1, 
Y{L) — —1/2 and Y{Hd) — —1/2 and thus, the superpotential term is invariant under LV(1) transformations. How- 
ever a non-zero vacuum-expectation- value of the scalar field corresponds to a electric charge breaking minimum. 
In the same sense also color breaking minima arise from the potential in a completely analogous way. The undesired 
global minima of the scalar fields can be translated into charge- and color-breaking bounds (CCB). Let us sketch the 
bounds found on the A-parameter, where we follow closely Ref. |106j . We start with a generic trilinear superpotential 
term W = A(^i02 03 with a corresponding scalar part W = X(j)i(j>24'3- As shown in App. [c]we get from this superpo- 
tential term a physical potential by collecting the F-terms, D-terms as well as the corresponding soft supersymmetry 
breaking terms, yielding mass terms and trilinear A-parameter terms. The physical potential thus reads 



^0(01,02,03) =iAp (i0in02p + i02n03i'+ i0in<^3 

+ ffa iYl\c^i\' + Yl\c^2\' + Yl\cj^,\'] 



102 



'03 



103 



(5.3) 



(AA010203 + c.c.) . 



Here we denote with Y^,, i = 1,2,3 the eigenvalues of the gauge group generators with adjoint index a and corre- 
sponding couplings ga, originating from the D-terms. Since the initial superpotential term (along with its derived 
Lagrangian terms) is supposed to be gauge invariant we have to have Y^^ — —{Y^_^ + Y^^). Now we are looking for 
the global minimum of the potential V. To this end we examine the directions in field space with = 0i = 02 = 03- 
In this direction in field space, the so-called D-flat direction, the quartic D-terms of the potential vanish, and do 
not protect the potential from a stationary solution for non-vanishing fields. Moreover, in the D-flat direction the 
potential becomes very simple, that is 



y^(0) = 3|A|2|0|4 + (ml^ + + J |0|2 _ 2A\<t>^ 



(5.4) 



Of course we have the desired stationary solution for = with V0(O) = 0. In order to avoid a vacuum for non- 
vanishing fields and thus to avoid a charge breaking minimum, there has not to be a stationary solution with a negative 
potential value. This eventually restricts the A-parameter not to be too large, that is 



^2 < 3 (r 



'03 J 



(5.5) 



A further approximative constraint arises from the Higgs potential with respect to the Higgs singlet [501 11021 1106] . 
The dominant singlet dependent part of the Higgs potential (3.8 1 reads 



Vs{S) = n'S^ 



(5.6) 
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where the ellipsis denote terms which have a lower dependence on S. Generally, we want to have a non-vanishing 
/i-term, with /i = Xvg, requiring a non- vanishing vacuum-expectation- value Vs- That is, in this case we have to have 
a minimum with a lower potential value compared to the symmetric minimum at Vs{0) ~ 0. This immediately 
translates into the approximate parameter condition 



ml<lAl. (5.7) 



Note, that this relation is only approximatively valid, since terms in the potential (5.6) are neglected 



The authors Miller, Nevzorov and Zerwas [104] study the impact of the vacuum stability constraints as well as the 
experimental constraints on the Higgs-boson sector parameter space. They derive simple analytic expressions for the 
physical Higgs masses in the CP-conserving case, taking the one-loop contributions to the Higgs-boson masses from top 
and stop loops into account. Using a power expansion in both l/tan(/3) and 1/Ma, where Ma is the upper left entry 



in the pseudoscalar mass squared matrix (3.271, they find analytic expressions, approximative valid for moderate or 
large values of tan(/3) and a larger scale Ma- Three distinct regions are considered with respect to the PQ-symmetry 
breaking parameter n. The region with vanishing k corresponds to the PQ-symmetric NMSSM, the region with large 
values of k, denoted as the NMSSM with strongly broken PQ-symmetry and the region with small values of k, that is 
K ^ 1 , denoted as slightly broken PQ-symmetry. Of course, the negative searches for a light axion on the one hand 
and the requirement of absence of a Landau-pole through renormalization group equations for k up to the GUT scale 



favor the slightly broken PQ-symmetry scenario (see Sect. 5.1 1. As an example, the Higgs-boson masses are plotted 
in a slightly broken PQ-symmetry scenario, namely with k = 0.1 in Fig. |8] The other parameters choices are given 
in the figure caption. Also the strong experimental constraint on the niA parameter is indicated in this figure. Note, 
that niA in general is not a CP-odd Higgs boson mass in the NMSSM, but the upper left entry in pseudoscalar mixing 



matrix squared; see (3.281. The authors point out, that the spectrum, based on their assumptions, in the NMSSM 




FIG. 8: The one-loop Higgs-boson mass spectrum as a function of Ma, the upper left entry m the CP-odd mixing matrix. The 
remaining parameters are fixed to n — 0.1, depicted as slightly broken PQ-symmetry, Vs = 3v, tan(/3) = 3 and = —100 GeV. 
The arrows denote the region allowed by LEP searches with 95% confidence. Figure taken from llO-j'j - 

is quite different from what is to expect from the MSSM. So, even if some of the Higgs bosons are too heavy to be 
detected, the discovery of the lighter Higgs-bosons may allow to distinguish the NMSSM from the MSSM. 

5.2. Experimental constraints 

1. Collider constraints 



The collider constraints give very clear bounds originating from very different observations like the Z-boson width as 
well as the direct searches of supersymmetric partner particles in e+e~ collisions at LEP. Here we discuss the bounds 
which are applied in the NMHDECAY program |107l 1108] : see App. |B]for an overview of some currently available 
program tools with respect to the NMSSM. 
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A strong constraint arises form the Z-boson precision measurements at LEP |109j . A possible contribution for the 
decay into the Ughtest and stable neutralino pairs Z —> XiXi with mass rn^o < mz/'2 may spoil the measurement of 
the invisible decay width r™'''^'*' = T'/ From a resonance scan the total width r^"' = 2.4952 ±0.0023 GeV 

(along with the Z-boson mass, mz — 91.1876 ± 0.0021 GeV) is known. The visible decay width consists of decays 
into charged leptons (r'^ '~ = 83.984 ± 0.086 MeV, with / = e,n,T) and hadrons {T^""^ = 1744.4 ± 2.0 MeV). From 
this the invisible decay width is measured as Y^^^^^^^^° — 499.0 ± 1.5 MeV. This is about what we expect from the 
SM, which predicts an invisible part, consisting solely of neutrinos as p™'''^''^' _ 501.3 ± 0.6 MeV. This means a 
possible contribution of decays into neutralinos should not exceed the measurement to much, 

r(^^X?X?) <2MeV. (5.8) 

Neutralinos may be produced in pairs at LEP via s-channcl Z-boson or Z-boson 7 interference or via t-channel exchange 
of a selectron. A further constraint arises from searches for pair production of non-LSP neutralinos with subsequent 
decay. These processes were searched for at LEP up to center-of-mass energies of 208 GeV [110 . From this negative 
search results the limits on the cross sections 



(7(e+e- ^ X?X°) < 0.01 pb , 
aie+e- ^X°X°) <0-lph 



r,o,m , (5-9) 



with i,j — 2,..., 5 can be deduced, viable for a sum of neutralino masses in the final state not exceeding the 
center-of-mass energy of the colliding electrons. 



Charginos could have been produced in pairs at LEP via s-channel exchange or via t-channel exchange of a sneutrino. 
The negative search can be translated into a minimal chargino mass limit of |lllj 

m^± > 103.5 GeV. (5.10) 

The production of charged Higgs bosons in pairs e+e^ H^H^ was investigated at LEP |112j for general SM 
extension with two Higgs doublets, like the MSSM or the NMSSM. Assuming main decay channels of the charged 
Higgs bosons cs and iJ+ t'^v^. and charged conjugated decay products for the H~ , a combination of all 

four LEP experiments ALEPH, DELPHI, L3 and OPAL with center-of-mass energies up to 209 GeV, give a lower 
mass limit of 

m|> 78.6 GeV. (5.11) 

There are constraints from the negative searches of the lightest CP-even Higgs boson Hi based on CP-conserving 
NMSSM Higgs sector studies. The dominant production at LEP proceeds via Higgs-strahlung off a s-channel Z-boson 
(e+e" Z* ^ ZH). This Higgs-boson production channel was searched for in various subsequent decay channels. 
The bounds derived in this way depend on the Higgs-boson mass. Here we mention the subsequent decays H —t bh 
and H t+t~ [75], H jj, with j denoting a jet |113l 1114] , ^ 77 |115j . H invisible, that is into LSP 
neutralinos |1161 1117j . H ^ X, that is independent of the decay product X in e+e~ ZX. The latter channel is 
accessible by studying the recoil mass spectrum in Z —f e^e~ and Z —> fJ.'^ fJ.~ events and by searching for X — s- e+e~ 
or X ^ 77 and Z ^ lyp ^ I116j . H — > AA with subsequent decay AA 4j, AA 2jcc, AA 2jr+r~, 
AA T^T^T^T^ , AA cccc, AA —f t^t^cc |118j . Also it was looked for the associated production of two Higgs 
bosons, e+e" —f HA, with subsequent decay HA bbbb, HA ^ t+t~t+t~, and HA AAA 3(66) |119j . 



28 



2. Muon anomalous magnetic moment 



The muon anomalous magnetic moment is a quantum effect which is very valuable in studying the possible deviations 
from the SM. For a recent review we refer the reader to Stockinger |120j . 
The magnetic moment of an electron or muon is 



9\ w- IS 

2m 



(5.12) 



with TO, e and s mass, charge and spin of the electron or muon and g the gyromagnetic ratio, which is 5 = 2 following 
from the Dirac equation at the classical level. Quantum effects lead to a deviation of this value, called the anomalous 
magnetic moment 



1 



(5-2) 



(5.13) 



The one-loop QED correction was first calculated by Schwinger and is a = ^ |121j . The contributions to the muon 
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FIG. 9: Contributions to the muon magnetic moment. From left to right the first diagram gives the contribution to the 
normal gyromagnetic ratio g = 2. The deviations to g = 2, that is the anomalous magnetic moment come from the quantum 
corrections shown in the five right diagrams. The second diagram gives the leading QED Schwinger contributing ^ to a in 
(5.131, the third and forth diagrams are the leading electroweak contributions. The fifth diagram shows the leading hadronic 



contribution via vacuum polarization and the last one the hadronic light-by- light contribution; see the review |120) . 

magnetic moment are shown in Fig[9j Generally, loop contribution from heavy particles with mass M are suppressed 
by a factor rn? /M"^ |120| . This makes clear that the anomalous magnetic moment of the muon, a^, is enhanced by a 
factor (to^/toe)^ k, 40,000 compared to the electron with respect to these contributions. 

For muon anomalous magnetic moment calculation in the SM we refer to the more recent publications mentioned 
by the Particle Data Group, |122l 11231 11241 WI^ 11261 [T?7] . For the corresponding calculation for the muon anomalous 
magnetic moment in SUSY models see [UHl HIS 1131 USD HSll ESS]. 

Let us briefly present the current results. The Brookhaven g — 2 experiment E821 |134[I135] measures an anomalous 
magnetic moment of the muon of 



a^^^P = (116, 592, 080 ± 63) x 10 
compared to the prediction of the SM, see the review 
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al^^ = (116, 591, 785 ± 61) x 10" 



Thus we get a deviation of 



5a,, 



SM 



(295 ± : 



X 10" 



(5.14) 
(5.15) 
(5.16) 



where we added the errors quadratically. If we trust the measurement as well as the SM prediction, this is a 
more than 3-(t deviation. The situation is presented by the Particle Data Group [109] and displayed in Fig. 10 
In this figure the deviations of the predictions from the BNL measurement are shown. Note that the prediction 
based on the r data which enter the hadronic vacuum polarization contribution do not deviate very much from the 
measurement, whereas all other predictions based on e^e~ data show a large deviation. With an E821 upgrade 
proposal, called E969 |137j , the error is aimed to be at least halved and with the same improvement of the theoretical 
predictions |136j . the 5-cr discovery limit may be reached. 

Of course supersymmetry changes the predictions of a^, since we get also quantum correction contributions from the 
superpartners. The outcome depends on the masses and couplings, not to mention the supersymmetric model under 
consideration itself. In any case, the experimental result (5.141 gives a strong constraint on the model. Additional 
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DEHZ (e*e"-based) 
-275 + 56 

DEHZ (T-based) 
-80 ±63 

HMNT (e*e"-based) 

-276 ±51 

J (e*e"-based) 

-287 ±68 

TY (e*e"-based) 

-274 ±59 



BNL-E821 (average) 

0±63 
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FIG. 10: Recently published predictions of the SM anomalous magnetic moment of the muon a^, subtracted from the mean 
experimental value from Brookhaven E821, as presented by the PDG |lU9j . The dots are predictions based on e^e~ data, the 
small triangle the prediction based on t data and the experimental measurement is given by a little square. The predictions 
are taken from ^22i ll23| il24j il25l [1261 [122] 




FIG. 11: Additional Feynman diagrams contributing to the anomalous magnetic moment of the muon in the MSSM and 
the NMSSM. 



Feynman diagrams to lowest order contributing to the anomalous magnetic moment in the MSSM and the NMSSM 
are shown in Fig. 1 1 For large tan(/3) the dominant supersymmetric contribution comes from the chargino-sneutrino 
loop diagram and is approximately given by (138j 



SUSY I 



a 



tan(/3) 



(5.17) 



with a the fine structure constant and M the heavier of the masses in the loop, that is of chargino and sneutrino. 

The additional contribution to in the constrained NMSSM (cNMSSM) (The cNMSSM is introduced in App. [Ef 
is shown in Fig. 
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in dependence on the unified gaugino mass where the unified scalar mass mo is fixed to 

zero |1|12J. We see in this plot that the cNMSSM gives for M1/2 % 1 TeV a contribution in agreement with the 
observation at BNL. 



It is worthwhile to add some critical remarks: first of all the experimental value with this high precision of (0.54 
ppm) is not confirmed by an alternative laboratory. Secondly, the uncertainty of the predictions is obvious from 
the deviating results based on the one hand on e^e~ data and on the other hand on t data. Moreover, even if the 
theoretical predictions of the hadronic leading contribution may be under control, the hadronic light by light scattering 
contribution (HLLS) is more or less estimated with a rather large contribution of a}^^^^ = 110 ±40 |136LI139] . Some 
effort is done to make this prediction more reliable; see |120j . 
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FIG. 12: (5a^^^*' m the cNMSSM as a function of the unified gaugino mass Mi/2- The darker and brighter shaded areas show 
the 1-a respectively 2-a deviations from 5a^ — ajf — afi_'^' . Figure taken from \10Slj . 



3. B-meson decay 



The flavor changing neutral current decay T(b — > 57) is very sensitive to physics beyond the SM. The reason is, that 
this decay can only proceed via loops in the SM, thus is suppressed. Hypothetically postulated new particles may 



contribute to the decay. In supersymmetry there are lots of contributions, two examples are shown in Fig. 13 that 



is a neutralino-sbino loop and a Higgs-boson-top loop. However a quite good agreement of experimental observation 

W X° H- 

: t t \ s 

\l \l \l 

FIG. 13: Example diagrams contributing to the b — > 57 decay. The left diagram is a SM contribution, where also the CKM- 
matrix entries Vth and V^* are depicted. The mid- and right-diagrams show contributions occurring in supersymmetry via a 
neutralino-sbottom loop and a Higgs-boson-top loop, respectively. 



and SM prediction of the & — > 57 decay is found not allowing for a large additional total contribution. But, on the 
other hand, stringent constraints are put on every new model which postulates new particles which may contribute 
to this decay channel. 

The experimental world average for the branching ratio (for a minimal photon energy of i?min = 1-6 GeV) is 
presented by the Heavy Flavor Averaging Group |140] based on the measurements of BABAR |141j . Belle |142j 
CLEO [143] : 

BR°''P(& ^ S7) = (3.52 ± 0.25) x 10^"* . (5.18) 
In contrast, the corresponding SM prediction |144j is, including the charm loop |145j as well as two- loop contributions 

BR^'^(6 S7) = (3.73 ± 0.30) x 10^'' . (5.19) 

That is, within the errors we have good agreement of the experimental data with the SM prediction. Thus, additional 
contributions from extensions of the SM are very restricted in case they contribute to the b sj decay. 
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5.3. Cosmological constraints 



Supersymmetric extensions of the SM predict a stable lightest supersymmetric particle (LSP), as long as matter 
parity (or R-parity) is conserved. This LSP may account for the observed cold dark matter in the Universe. In 
the early Universe the LSP would be present in large numbers in thermal equilibrium. As the Universe has cooled 
the LSP's only may reduce their density by pair annihilation and co-annihilation |157l 11581 1159] . As their density 
decreases, the chance to find another particle to annihilate decreases and at a certain point the comoving density 
can become constant ("freeze out"). The remaining relic density can be calculated depending on the considered 
model. In case the relic density of the LSP is identified with the cold dark matter candidate, a comparison with the 
cold dark matter observations constrain the model. We distinguish between direct detection of the relic density LSP 
particle in a detector and the indirect detection via astrophysical observations. 



1. Indirect dark matter detection 



The most stringent bound comes from the cosmic microwave background (CMB) observation performed by WMAP. 
The observed CMB is compared to the predictions of the standard model of cosmology, the so-called A-cold-dark- 
matter model (ACDM). This model depends on six parameters, the Hubble parameter Hq, the baryon density ^IbfiQ, 
the cold dark matter density ^Ic^q, the dark energy density JIa, the scalar spectral index rig, the optical depth r, as 
well as the variation of the spectral index and the curvature perturbation A^. By adjusting the parameters of the 
ACDM, the prediction of the CMB power-spectrum can be fitted to the observation. Taking into account the recent 
five year WMAP data, combined with measurements of type la supernovae and baryon acoustic oscillations in the 
distribution of galaxies, the current 1-cr limits on the ACDM parameters are given as |160j 

^bhl = 0.02267l:o;ooo5g , 
f^c/io = 0-1131 ± 0.0034, 
r^A = 0.726 ±0.015 , 

5.20 

= 0.960 ±0.013 , ^ ' 

T = 0.084 ±0.016 , 
= (2.445 ± 0.096) x 10"^ at fc = 0.002 Mpc"^ , 

where we are in particular interested on the cold dark matter density fic^o constraint. Note that the Hubble constant 
is defined in terms of ho, that is Hq = 100 ho km/s/Mpc. 

Now let us briefly sketch the calculation of the cold dark matter density in supersymmetric models. This relic 
density of the LSP as cold dark matter candidate, supposed to be a neutralino, can be computed by the solution of 
the continuity equation [1581 1159] 

f^{n<,R') = -{a^nnv) (n| - (nl^f) . (5.21) 

This continuity equation states that the number of neutralinos, with density n^^, in a comoving volume, is governed 
by the competition between annihilations and creations. The creations depend on the density of neutralinos in 
equilibrium n-'^. It is assumed that the annihilating particles have non-relativistic relative velocity v. Further, it is 
assumed that the annihilating particles as well as the annihilation products are maintained in kinetic equilibrium with 
the background thermal plasma through rapid scattering processes |161l I162j . The model-dependent input is given by 
the thermally averaged annihilation cross section ((Tannt^) = a ± bv^ with parameters a and b. With this annihilation 
cross section the cosmic neutralino density can be deduced |163l 11641 1165] 

2 imxlO^Xfr 1 



with Alp the Planck mass and 5, ~ 81 the effective number of degrees of freedom at the freeze-out temperature 
Tfr and Xfr — vny^jTjr-, where is the neutralino mass. The annihilation cross section in models like the MSSM 
and the NMSSM can be computed from s-channel Higgs and Z-boson exchange into fermions and gauge bosons, 
and from heavier neutralino or chargino t-channel exchange with gauge bosons in the final state. In particular 
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for s-channel resonances, rapid annihilation can take place, reducing the relic density strongly. In the NMSSM, 
in contrast to the MSSM we have additional Higgs-boson resonances, contributing to the annihilation. In the 
case of light Higgs bosons in the NMSSM also new decay channels of neutralinos decaying into a Z boson and 
a Higgs boson or two light Higgs bosons are available jl66l 11671 11681 11691 11701 I171j . Contrary, a large singlino 
component of the LSP reduces the annihilation cross section, since singlinos do not couple to gauge bosons at tree level. 



In the investigation of Belanger et al. jlTlj the parameters of the NMSSM are scanned in particular with view on 
the dark matter constraint in (5.201. In this publication the program package micrOMEGAs |172j is employed, which, 
starting from the particle spectrum and couplings, computes the relevant annihilation and co-annihilation neutralino 
cross sections and eventually the relic density. The particle and coupling spectrum is provided by the NMHDECAY 
program |108j . Assuming gaugino mass unification at the GUT scale (Mx/2), the independent parameters are A, k, 
tan(/3), /i = \Va, A\, A^. The authors find parameter space corresponding to a large singlino component of the 
neutralino LSP. This parameter space passes also the current collider constraints. Some representative examples are 
given in Tab. Ill [172 . In the upper part the initial NMSSM parameters are given, followed in the mid part by the 



parameters 


1 


2 


3 


4 


5 


6 


A 


0.6 


0.24 


0.4 


0.23 


0.31 


0.0348 


K 


0.12 


0.096 


0.028 


0.0037 


0.006 


0.0124 


tan(/?) 


2 


5 


3 


3.1 


2.7 


5 


\vs [GeV] 


265 


200 


180 


215 


210 


285 


Ax [GeV] 


550 


690 


580 


725 


600 


50 


A^ [GeV] 


-40 


-10 


-60 


-24 


-6 


-150 


M2 [GeV] 


1000 


690 


660 


200 


540 


470 



spectrum 



m^o [GeV] 


122 


148 


35 


10 


15 


203 


t/f 3 + VIa 


0.12 


0.29 


0.12 


0.03 


0.06 


0.02 




0.88 


0.69 


0.87 


0.95 


0.94 


0.96 


m^o [GeV] 


259 


199 


169 


87 


182 


214 


m^l [GeV] 


258 


193 


171 


139 


196 


266 


mn' [GeV] 


117 


116 


36 


22 


34 


115 


^13 


0.88 


0.04 


0.98 


1.00 


1.00 


0.04 


[GeV] 


128 


158 


117 


114 


113 


163 


-R23 


0.11 


0.96 


0.01 


0.00 


0.00 


0.96 


ruA, [GeV] 


114 


59 


56 


18 


15 


214 


-R45 


0.99 


1.00 


0.99 


1.00 


0.99 


1.00 


dark matter 


study 














0.1092 


0.1179 


0.1155 


0.1068 


0.1124 


0.1023 


channels 


ha (73%) 
VV (13%) 
ZH (8%) 
HH (3%) 
m (2%) 
11 (1%) 


VV (75%) 
ha (17%) 
HH (5%) 
ZH (2%) 


qq (65%) 
U (35%) 


qq (93%) 
// (7%) 


aa (92%) 
qq (7%) 
11 (1%) 


X2X2 ^ X (81%) 
X?X2 - X (15%) 
X?xf - X (2%) 
qq (2%) 



TABLE III: Parameters with a LSP neutralino with a large singlino component {Ut-s) respecting the current collider as well as 
the cosmological dark matter constraints as pres ented in I171j. The neutralino mass mixing matrix U is defined in (D27I, the 
Higgs boson mass mixing matrix R is defined in (3.221. Here a CP-conserving Higgs sector is assumed. 



particle spectrum and mixings and in the lower part by the dark matter study. The annihilation and co- annihilation 
channels are given which lead to the required relic density ^^ali^ in the 3-(T range of fic^O' ^ given in (5.201. For 



the parameter points 3-5 there are very light LSP neutralinos which have, following from the neutralino mass mixing 
parameter C/^g (the matrix U is defined in (D27l), a large singlino component. In this case the dominating annihilation 
proceeds via Z- and Higgs-boson resonances. The parameter points 4 and 5 correspond to a light scalar Higgs boson 
which is dominantly a Higgs-boson singlet, as can be seen from the Higgs-boson mass mixing parameter « 1 
(the matrix R is defined in (3.22 1). In this case the lightest non-singlet SM-like Higgs boson has to pass the LEP 
bound. This in turn drives tan(/3) close to 3, where the mass mjj^ is maximized. In Fig. 14 the impact of the WMAP 



constraint on the parameter space is shown. The contours for different values of f2^o/iQ are given, where the black 
region corresponds to the WMAP constraint, that is 0.0945 < fi-j^o/iQ < 0.1287. The chosen independent parameters 
are described in the figure caption. For further details see |171j . 
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H [GeV] 



FIG. 14: Contours for Q^ohg — 0.02, 0.0945, 0.1287, 1 in the fi-M2 plane, as presented by \1 71f . The remaining parameters 
are fixed to X = 0.1, k = 0.1, tan(/3) = 5, Ax — 500 GeV, ~ 0. The gaugmo masses are assumed to umfy at the GUT scale, 
corresponding at the electroweak scale to M2 — 2A/i and Mg — 3.3M2. The hatched region is parameter space which is excluded 
by LEP due to searches for chargino pairs (5.101; see 1171^ for further details. 



2. Direct dark matter detection 



If the dark matter consists of unidentified particles, the earth should be passing through a flux of these particles 
which constitute the dark halo of our Milky Way. Direct detection of dark matter may then proceed via elastic scat- 
tering from a target nucleus. The nuclear recoil could then be detected in an appropriate detector. The corresponding 
event rate depends on the dark matter density in the solar vicinity and the model-dependent neutralino-nucleus elastic 
scattering cross section, under the assumption that neutralinos make up dark matter. For a fixed dark matter density 
(which is, following cosmological arguments about pdm ~ 300 MeV/cm'^) the flux of dark matter particles through a 
detector is obviously inverse proportional to its masses. 

For investigations of direct neutralino dark matter detection in context with the NMSSM we mention the approaches 
[1701 11731 11741 11751 I176j . On the experimental side we draw the attention to the experiments of the DAMA col- 
laboration [T77], CDMS [T75], EDELWEISS [T75], ZEPLIN I [THO], as well as the XENONIO experiment [HI]. Up 
to now there are no confirmed detections of any weakly interacting massive particle (WIMP). The sensitivities of 
some direct detection experiments are shown in Fig. 15 depending on the mass of the neutralino |176j . The dots give 
predictions of the NMSSM varied in certain parameter ranges, as indicated in the figure caption. These points pass 
the LEP/TEVATRON, muon anomalous magnetic moment and BR{b 57) constraints as discussed above in this 
section. The grey dots have passed the relic density constraint 0.1 < ilhg < 0.3 and the black dots the much stronger 
constraint from the three-year WMAP data 0.095 < HHq < 0.112. Obviously there is parameter space left which 
passes all the experimental and indirect/direct cosmological constrains in the NMSSM. As indicated in the figure, 
future direct dark matter detection experiments may give strong constraints of viable parameter space. 



3. Baryogenesis 



The absence of antimatter, that is the observed baryon asymmetry in our Universe may be generated by strong 
electroweak phase transitions (EWPT) of first order EH ESI |27] . Alternatives to EWPT like GUT baryogene- 
sis |182j . baryogenesis via leptogenesis |183j and the Affleck-Dine mechanism |184j are not discussed here. For an 
overview, see for instance |185j . One of the attractive properties of EWPT is that this mechanism is understandable 
in terms of physics at the electroweak scale, accessible at current colliders. Our focus lies on the discussion of EWPT 
with respect to the NMSSM. We start with a brief sketch of EWPT and show how parameter constraints are deduced 
in the SM, the MSSM and in the NMSSM. 

Phase transitions are generally characterized by an order parameter which suddenly changes at a critical temperature 
Tc- In case the order parameter has a discontinuity at Tc, the transition is called a first-order phase transition. In 
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FIG. 15: Scatter plot of the neutralino-nucleon cross section as a function of the neutralino mass, here denoted by for 
tan(/3) = 5, and the remaining parameters varied in the ranges 0.01 < A, k < 0.7, 110 GeV < M2 < 430 GeV, —300 GeV < 
< 300 GeV, -800 GeV < < 800 GeV, 110 GeV < VsX < 300 GeV. All points pass the LEP/TEVATRON collider, 
muon anomalous magnetic moment and BR{b — > 57) constraints and have a relic density in the range 0.1 < flchg < 0.3 (grey 
dots) or 0.095 < flchg < 0.112 (black dots) reflecting the three year WMAP data. The sensitivities of experiments are given by 
full lines, whereas the dashed lines give the sensitivities expected in the future. Figure taken from I170j . 



contrast, if the order parameter changes continuously at T^, it is called a second-order phase transition. In the 
electroweak phase transition the order parameter is identified with the Higgs vacuum-expectation- value. Assuming a 
first order electroweak phase transition we have with decreasing temperature a transition from the symmetric phase 
with a global minimum at the origin for T > Tc to a symmetry breaking phase with a global minimum with a 
vacuum-expectation- value at w ^ for T < Tc. For T — Tc we have degenerate minima at w = and v — Vc. 
It has been shown that a strong first order phase transition is needed [186} 1187) , that is 

g>^«l (5.23) 

in order to avoid a wash-out of the generated baryon asymmetry after it is generated. 

Following closely the argumentation found in [ 1881 1189j let us consider a generic Higgs potential 

V' = m^<tP - act)^ + (3(t)^ (5.24) 

with a generic scalar Higgs field 0. From this potential we see that we get a first order phase transition, that is 
degenerate minima with a symmetry breaking minimum at Wc, for 

(5.25) 

Obviously we need a non- vanishing parameter a, that is a cubic term in the potential in order to get the wanted 
symmetry breaking minimum. 

In the SM and in the MSSM we do not have a cubic term at all. Nevertheless, an effective cubic term is generated 
by thermal loops [TSOl EI] 

Kff =V' + y^oop (5.26) 

with 

^Soop = rET/(g,ln(l± (5.27) 



35 



with i running over all fermions (upper sign) and bosons (lower sign), which couple to the Higgs boson. By mi{4>) are 
the corresponding field-dependent masses denoted, that is the masses depending on the Higgs field, before electroweak 
symmetry is broken. For high temperature T the effective potential can be expanded as [188 



-m,((/)) In 



m1{4>) 



-1 
4 



O 



T 



where the upper line is valid for bosons and the lower one for fermions. The constants are 

'|+21n(47r)-S 
Cboson - 21n(4) 



_ j I + 2 ln(47r) - 2'yE , for bosons 



for fermions 



(5.28) 



(5.29) 



Obviously, we get for each boson in (5.28) an effective cubic term. For instance, for the field-dependent charged 
M^^-boson masses mY/±{4') = g24>/V^ a cubic a-term in (5.241 is generated and first order phase transitions become 
possible. 



In the SM however, these loop effects are too small, thus give only a too small a-term. From the condition (5.231 
one can deduce an upper bound on the Higgs-boson mass and gets |192j 



m 



If < 32 GeV, 



which is excluded from the lower bound on the Higgs-boson mass from the LEP experiment, 
95% CL [75 . The conclusion is, that this mechanism of strong first order EWPT is not possible in the SM. 



SM, cxp 



(5.30) 



> 114 GeV at 



In the MSSM the situation is quite different. Let us sketch briefly the situation in this model. In Sect. 3.7 we have 
already seen that the tree-level mass of the lightest Higgs boson, (m^^^^)^ < m| cos^(2/3), is ruled out experimentally 

by LEP from the lower bound of rri^^^^' > 92 GeV (H]. However, as we have discussed, there are quite large 
quantum corrections to this tree level mass. The largest contribution to the quantum corrections arise from top and 
stop loops, which couple strongly to the Higgs-bosons. We recall the one- loop radiative contribution to the lightest 
Higgs-boson mass (3.37), 



(A? 



MSSMn2 



In 



mi m,i 



(5.31) 



This correction allows for larger Higgs-boson masses passing the current LEP limits. From the logarithmic term 
we see that large radiative corrections require at least one of the masses m^^, m^^ to be very large. From precision 
electroweak constraints we know that has to be the heavy scalar. This follows from the fact that couples stronger 
to the electroweak gauge bosons than t^, and a suppression of this loop contributions requires that is the heavier 
scalar. However, from the logarithmic dependence on m^^ of the mass shift (5.31 1 it follows that a substantial shift 
of the Higgs-boson mass requires a very large mass m^^ . 



Now we come to the thermic loops which potentially give a cubic term in (5.28 1 in the potential required for a strong 
first order EWPT [Ml [TMJ [TM CSi] . The field-dependent stop mass m^^ 
the stop mass itself into account. 



reads, taking also thermic corrections of 



TO? 



mi + \yA'\Hy + cnT 



(5.32) 



with CR = 4/9g| + l/6|yt|2(l + sin^{l3)). A strong EWPT requires a large cubic term proportional to and this 
in turn requires to| ~ —crT^, that is a negative soft-breaking parameter to~. For vanishing temperature this means 
that we have to have a light i^j-mass below the top-quark mass. Altogether we find a mass hierarchy in the MSSM 
TOj^ < TOt << TO(^ , which is considered to be unnatural, since it requires fine-tuning. In Fig. 16 the experimental 
95% exclusion regions in the stop-neutralino mass plane are shown as presented by |196j . We see that for a light 



neutralinos to^ < 45 GeV stop masses below 120 GeV are excluded. Thus, with respect to (5.31) we certainly need 



corrections of the order of 30 GeV in the MSSM, that is rather large heavy scalar top mass TOj^ compared to to^^, 
that is large fine-tuning. 



In the NMSSM the strong electroweak phase transition requirement (5.23) is generically much easier fulfilled since 



already at tree level we have cubic terms in form of additional soft-breaking terms |197j . 

Kioft, trilinear = \Ax{HleHd)S + A^S^ + C.C. 



(5.33) 
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rrij (GeV/c^) 



FIG. 16: 95% exclusion regions in the stop-neutralino mass plane for the MSSM. Figure taken from fl96f . 

which are absent in the SM and the MSSM. Note that a cubic term in the generic Higgs potential can not be 
constructed in a gauge invariant way from Higgs doublets, since gauge invariance requires the Higgs doublets to 
occur in a bilinear or quartic form. Thus, an additional singlet Higgs field S, as introduced in the NMSSM, opens 
the possibility of cubic terms without relying on quantum corrections. Moreover, the tree-level mass limit on the 
lightest Higgs-boson mass in no longer valid in the NMSSM. 



There was some effort spent in systematical numerical parameter scans in order to find parameter space fulfilling 
the strong EWPT constraint (see for ins tance ^TM fTWl [TM [TOS] ) . Let us start with the study of Pietroni ^7\. The 
Coleman- Weinberg one-loop corrections (3.35 1 as well as the non- vanishing temperature contributions to the potential 
were taken into account |1901I191| . The following assumptions were imposed. 

• Common soft-breaking masses ttiq = = Af| = 1 TeV. 

• Fixed values for X'^{mz) = 0.274, yt(mz) = 0.97 and = l/2X'^. 

• A minimal value for the singlet VEV Xv^ < 45 GeV, following approximatively from the experimental lower 
bound on the chargino mass. 



The ratio tan(/3) is fixed to the values 2 and 10 and for the Higgs-boson mass parameters it is assumed that 



mjj > 0, < and m| > 



2 - n ^2 ^ ^ ^ ^^^^ ^ ^ 

NMSSM 



From the parameter scan it is found that the allowed parameter space leads to an upper bound of < 170 GeV 

for the lightest Higgs boson, far beyond the LEP limits. 



In the work of Bastero-Gil et al. [189J also the parameter space is examined targeting on the requirement of strong 
first order EWPT. With different assumptions for the soft-breaking parameters as well as for the experimental LEP 
bounds it is shown that there is parameter space available in the NMSSM which also respects the LEP Higgs-boson 
mass bound. 



Let us also mention the paper of Davies et al. [TijS ] . The authors report that in a parameter scan passing appropriate 
constraints about one half of parameter points pass the strong first order EWPT condition. Under certain assumptions, 
not given here, they find masses of the lightest Higgs boson not exceeding 120 GeV. In this study also a generalization 
of the superpotential W in (3.2 1 is introduced with an explicit /i-term and a linear term: 



W' = W + ^x{HjeHu) ~ rS . 
Of course, by setting the dimensional parameters /i and r to zero the original NMSSM is restored. 



(5.34) 



Eventually in the publication of Huber and Schmidt |199j emphasis is placed on a subtle point in context with 
first order EWPT. As they point out, for an EWPT it is required, that phase transition in Vd, and Vs have to 
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occur simultaneously. This means that all three vacuum-expectation- values should be of the same order of magnitude. 
If one has for instance a much larger value Vs compared to Vd and Vu, a cascade of phase transitions would occur 
with decreasing temperature. This would not correspond to a first order phase transition which needs a bump in 
the effective potential, requiring a trilinear term in the Higgs fields. The authors refer to the parameter scans of 
EUwanger et al. (see the discussion in Sect. 5.4 1, which predict a generic larger value for Vg compared to the doublet 
VEVs. In order to circumvent this problem it is proposed that the supplemented superpotential (5.341 should be 
taken into account. In this way a compatible vacuum-expectation-value Vs of the electroweak order arises. However, 
as is pointed out by the authors, the generic larger values of Vg arises in studies, where the assumption of unification 
of soft-breaking parameters at the GUT scale or other additional constraints are considered. As we will see in the 
next subsection, dropping the unification condition meets the criterion of a VEV Vs of the electroweak order also 
without a modification of the NMSSM superpotential. 
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5.4. Parameter scans 



A quite generic approach in order to constrain the parameter space is to scan over the parameter ranges of interest. 
For each parameter set in this scan the particle mass and couphng spectrum is derived. Then, the viabihty of 
this spectrum is checked, for instance with respect to the corresponding minimum structure: only a parameter set 
corresponding to a global minimum in the Higgs potential with the observed electroweak symmetry breaking is to be 
accepted in this scan. Moreover, the different experimental and cosmological constraints on the masses and couplings 
have to be passed. From the eventually gained allowed parameter space, which passes all constraints, restrictions 
can be read off. In this section we want to present only some of the various scans which have been performed in 
the NMSSM [SHJ \TM 12001 HHIl USD- Typically, they are based on computer tools which calculate the spectrum 
and apply various theoretical and experimental constraints. Due to the large number of parameters in the NMSSM 
typically not the full parameter space is scanned over for practical reasons. This is in particular true for the rich part 
of parameters originating from the soft supersymmetry breaking terms. Note that in the discussion of Higgs-boson 
phenomenology in Sect. |3.7| already some studies based on parameter scans were mentioned. Assumptions like the 
unification of parameters at a high scale like the GUT scale lead to a enormous reduction of available parameter space 
and thus simplifies the investigation. Therefore, most of the parameter scan studies are applied to the constrained 
NMSSM (cNMSSM); see Sect. |e] We start with briefly discussion some of these studies in the cNMSSM followed by 
a parameter scan in the general NMSSM without unification assumption. 

Let us first mention a series of investigations of the cNMSSM by EUwanger et al. [551 120011201] . For each of about 
one million points in this five-dimensional parameter space, the particle and coupling spectrum of the model at the 
electroweak scale is generated by the renormalization group equations. Then, the following constraints are applied: 



Large trilinear soft A-parameters may induce electrically charged or colored vacua; see (5.5). Since this is 
phenomenologically unacceptable, the corresponding parameter sets are discarded. 

The global minimum is searched for numerically and the global minimum is required to have non- vanishing Higgs 
VEVs Vd, Vu and Vg- In this, the effective one-loop Higgs potential is considered, taking the Coleman- Weinberg 



contributions (3.351 into account. The ratio of the VEVs Vd and is forced to comply with tan(/3) < 30. 



The top-quark mass, as derived from the particle mass spectrum, is forced to fulfill the experimental constraints, 
that is 168 GeV < mt < 192 GeV. Also it is checked that the negative search for charginos m^± is not violated; 



see (5.101 



It is found that low values of the parameters A and k are favored as well as a large value of the singlet VEV Vg- 
This means that the singlet is decoupled and the NMSSM mimics the MSSM under the unification assumption. 
Nevertheless, there remains parameter space where one neutralino has a large singlino component and thus substantial 
differences may arise compared to the MSSM. It is pointed out (see for instance [1021 ) that for small parameters 
A and k, the two parameters 171% and A^, change not very much between the electroweak and the GUT scale by 



the renormalization group equations, that is we get from (5.7 1 the unification parameter constraints TOq < 1/9Aq. 
Note that there is a subtle difference between the cMSSM and the cNMSSM: In the cMSSM small values of the mo 
parameter are disfavored since they lead to a charged slepton LSP, unacceptable as a dark matter candidate. In 
the cNMSSM contrary, a small parameter mo is favored due to a viable global minimum for non- vanishing Vs- The 
slepton LSP in the cNMSSM is avoided since the possibility of a singlino-like LSP arises in parameter space 



We would also like to mention the approach of Stephan [202| . This study is also performed in the constrained 
cNMSSM. The main difference to the previously discussed approaches is that in this work also the dark matter 
constraint, that is the relic abundance of the LSP, supposed to be a neutralino in the NMSSM, is taken into account. 
There are also some subtle points concerning some deviating theoretical and experimental constraints, which are 
applied. Starting with 5.5 • 10® points in the five-dimensional parameter space, about 4900 points pass the theoretical 
and experimental constraints but only 2000 points of these pass the cosmological constraint given by the relic neutralino 
abundance bound. Quite restrict explicit bounds are derived in this way and given already in the abstract of this 
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work: 



mH^< 140 GeV, 
m^o < 300 GeV , 
mg^ < 300 GeV , 
300 GeV < < 1900 GeV , 
200 GeV < rrii^ < 1500 GeV , 
350 GeV <mg < 2100 GeV . 



Scatter plots in the mg^-m^o plane are shown in Fig. 17 for the cMSSM as well as for the cNMSSM. Note that in the 
cNMSSM the selectron mass mg^ is stronger constrained than in the cMSSM. As mentioned above this arises from 
favored lower values of mo, whereas in the cMSSM small values of toq are disfavored. An upper bound of niQ « mg^ 
can be deduced based on the relic density constraint which requires to have nearly degenerate masses of the LSP and 
the NLSP |102j . From an approximation for the RGE running of the NLSP mass it is found that itiq < l/15M^y2; 



favoring low values of mo 



in the cNMSSM in contrast to the cMSSM. Also the strong effect of the application 



of the relic abundance dark matter constraint is demonstrated in this figure. 
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FIG. 17: Scatter plots in the constraint cMSSM (upper part) compared to the constrained cNMSSM (lower part) in the rUej^- 
m^o plane. The right part shows the scatter plots passing also the WMAP dark matter (denoted by d.m.) constraint. Figure 
taken from J20^ . 



In the work of Bastero-Gil et al. |189j . the MSSM is compared to the NMSSM. In this study emphasis is placed 
on the fine-tuning required in both models in order to comply with the theoretical and experimental constraints. 
Fine-tuning is quantified quite similar to (5.2), in the form 



A"^^'^ = max 



dm^ 



da,: 



(5.35) 



where denotes all the soft supersymmetry breaking parameters. In a parameter scan the unified scalar mass 
is fixed. Too = 100 GeV, as well as the A-parameters, ^^(to^) = 0, ^((Ag) — and the gaugino masses 
Mi(A(3) — M2(Ag) = 500 GeV, where Kq denotes the GUT scale. The gaugino mass Mg is varied in the range 
100 GeV < Mg{Ao) < 600 GeV, and the Higgs-boson mass parameter niHa is devoted to < mHai^a) < 1 TeV, 
with /i = Xvs < 0. In the NMSSM also the additional parameters X{Aq) = 1 and k{Aq) = 0.1 are fixed. Requiring 
further for the chargino masses to^± > 90 GeV, the values for the fine-tuning function A™'*'^ are given in Fig. 18 
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FIG. 18: Scatter plot of fine-tuning as defined in ( 5.35 1 as a function of the lightest physical Higgs-boson mass for tan(/3) = 3 
left and tan(/3) = 5 right. In this scan there are fixed values of mo = 100 GeV, At{A.G) ~ 0, Mi(Ag) = M2(Ag) = 500 GeV, 
— Xvs < and varied parameters Mg{A.G), mHji^c) and mH^i-Z^a) ■ The dark shaded regions shows the fine-tuning values 
corresponding to the NMSSM and the brighter shaded regions corresponding to the MSSM. The LEP limit on the SM Higgs-boson 
with mi^ < 108 mass is also shaded. Figure taken from I189f . 



for tan(/3) — 3 and tan(/?) = 5, respectively. Let us cite the authors: The plots are a striking demonstration that 
the physical Higgs boson can be heavier and involve less fine-tuning in the NMSSM compared to the MSSM at low 
values of tan f3. The large fine-tuning in the MSSM for low tan(/3) originates from the LEP bounds on the minimal 
Higgs-boson mass confronting the tree-level prediction of a mass below the Z-boson mass mz- Thus, large quantum 
corrections to the Higgs-boson mass in turn require a very large stop mass, that is large fine-tuning. On the other 
hand it is pointed out that low values of tan(/3) are favored in order to have a strong electroweak phase transition. 
The authors conclude that in the NMSSM the situation is much better with respect to the LEP bound on the 
minimal Higgs-boson mass, fine-tuning and baryogenesis via first order EWPT. 

In a recent study by Djouadi, EUwanger and Teixeira |102j, focusing also on the cNMSSM, very severe constraints 
are found. After applying the theoretical as well as current experimental constraints, including the new WMAP data 
on the relic density, there remains a quite narrow window of allowed parameter space. It is reported that in the 
cNMSSM, the gluino g turns out to be generically heavier than all squarks q. The NLSP is found to be generically 
stau-like and nearly degenerate with the LSP mass which is a neutralino, reflecting the constraints coming from the 
relic density. The supersymmetric partner particles eventually decay via the stau-like NLSP into the LSP neutralino. 
This opens the possibility of the observation of displaced vertices as discussed in Sect.|4] The particle spectra of the 
neutralinos, charginos, selectron and stau as well as the spectra of squarks and gluino are shown in Fig. [19) In this 
figures, the unified scalar mass is fixed to zero, mo — (also another choice of the unified scalar mass is discussed in 
the paper) . We see that the spectra of supersymmetric partner particles are highly restricted by the theoretical and 



experimental bounds, where also the current WMAP constraint for the relic density is applied; see (5.201. As the 
authors of this investigation have stressed, the measurement of one sparticle mass or mass difference would allow to 
predict quite accurately the complete sparticle spectrum in the cNMSSM. 

Let us also mention the approach of Bednyakov and Klapdor-Kleingrothaus [174J , targeting mainly at the direct 
detection of a LSP neutralino dark matter candidate. This direct detection refers to elastic scattering of a dark matter 
neutralino from a nucleus producing a nuclear recoil detected in an appropriate detector; see Sect. |5.3[ In this work 
no unification at the GUT scale of soft parameters is assumed, that is there is no need to apply the RGE's from the 
GUT scale starting point in this case. The scan is done over the eleven-dimensional parameter space 

Ml, Af2, tan(/3), Ws, A, k, mq^ , mg.^ , Au^^Ax, A^, (5.36) 

with rriQ-^ = toq2, the soft mass parameters with generation index. For simplicity reasons it is assumed that 

and niQ = m^, — — m^, = mg,. Moreover it is set A^ = Ad = A^ = 



2 2 9? 2 

rUn = ra- = m - = mi — m% aim 

W2 "1,2 di 2 1-2 '^i-^ V. 



for all generations except A^^. The soft supersymmetry breaking parameters are defined in (D16l. Negative results 



for the search for supersymmetric particles at LEP and at TEVATRON are included in the analysis. The resulting 
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FIG. 19: Spectrum of Higgs bosons tn the upper plot (here denoted with small letters), neutralino-, selectron- and stau-mass 
spectrum in the left plot and squark- and glumo-mass spectrum in the right plot. The unified scalar mass is set to tuq = 0. 



eneracy of the masses m-± = m^g and m.± — rrij^ 



: found. Plots taken from flU0j 



allowed parameter space is compared once without any cosmological constraint on the co-annihilation rate of the 
LSP and once taken the constraint 0.025 < ri-j^o/i§ < 1 into account. This cosmological constraint alone gives, as is 
reported, a 20% cut on the allowed parameter space. As is pointed out in their conclusion, under their assumptions, 
but without unification of soft supersymmetry breaking parameters at the GUT scale, there remain domains in 
parameter space where the lightest neutralino has a quite small mass, even as small as 3 GeV. In this work it is 
stressed that in a proposed germanium ^''Ge detector there are large event rates expected in the NMSSM. Let us 
remark, that with view on the current WMAP data, it would be very interesting to see the effect of the current much 



tighter bounds in this study; see (5.201 



Eventually let us remark that the results of parameter scans depend obviously strongly on the assumptions taken into 
account, like unifications of parameters. On the other hand are the results also strongly dependent on the constraints 
applied to the derived particle and coupling spectra. To summarize, we have seen that at least for the general NMSSM, 
without any unification assumptions, there is a lot of viable parameter space, even if the cosmological LSP cold dark 
matter constraint is taken into account. Nevertheless, even that the cNMSSM seems to be very restricted, it is a 
viable model with respect to all available theoretical and experimental bounds. 
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6. DETERMINING THE GLOBAL MINIMUM OF THE HIGGS POTENTIAL 



It is a non-trivial task to find the global minimum for Higgs potentials with a large number of Higgs fields. For 
instance in the NMSSM the Higgs sector consists of two complex electroweak doublets and one complex electroweak 
singlet, that is 8 real fields from the doublets plus 2 real fields from the singlet. The conventional approach based 
on the unitary gauge requires the global minimum to be found in a 7-dimensional field space. Typically, numerically 
methods are applied in such involved cases in order to determine the global minimum; in contrast, here we want to 
discuss an algebraic approach [63^. In this approach all stationary points of the tree- level Higgs potential are found 
and supposed, the potential is bounded from below, the global minimum is identified from the corresponding lowest 
value of the potential. The method is applied to the NMSSM, revealing a quite surprising structure of stationary 
points, that is minima, maxima, and saddle points with different behavior with respect to the symmetry breaking of 
the SU{2)l X U{1)y electroweak gauge group. 

The global minimum of the Higgs potential gives the expectation values of the Higgs fields at the stable vac- 
uum. Parameter values for the Higgs potential are thus considered to be acceptable only, if the global minimum of 
the Higgs potential occurs for Higgs field vacuum expectation values, which induce the spontaneous breakdown of 
SU{2)l X U{1)y to the electromagnetic U{l)em at the observed electroweak scale v « 246 GeV. 

Firstly, the tree level Higgs potential for general models with two Higgs doublets and an arbitrary number of 
additional Higgs singlets is considered. The first step is to notice that the potential is restricted by renormalizability 
and gauge invariance. Renormalizability requires at most quartic terms in the real fields in the potential. Electroweak 
gauge invariance restricts the possible doublet terms in the potential, since only gauge invariant scalar products of 
doublets can occur. Substituting the doublet fields by appropriate functions of these invariant terms, all gauge degrees 
of freedom can be eliminated from the potential and effectively reduce the occurring powers in the doublet terms. 
The method to base the analysis on quadratic gauge invariant functions was introduced already in context with the 
general two-Higgs doublet model |204l 1205] . 

The stationarity conditions form a non-linear, multivariate, inhomogeneous polynomial system of equations of third 
order. A systematic approach to solve these - in general quite involved - systems of polynomial equations is by 
a Groebner basis computation, well established in ideal theory [2061 12071 1208] . The Groebner basis was originally 
introduced to solve the ideal membership problem. Constructing this Groebner basis in an appropriate order of 
the monomials (the terms of the polynomials including coefficients), for instance the lexicographical ordering, and 
subsequent triangularization allows to solve the initial system of equations algorithmically for any finite number of 
complex solutions. The introduction of gauge invariant functions just avoids continuous gauge symmetries in the 
potential and the finiteness of the set of complex solutions can be easily checked within this algorithmic approach. 
Moreover, this approach guarantees that all stationary points are found. 

Then the method is applied to the NMSSM. For the computation of Groebner bases as well as the subsequent steps 
to solve the systems of equations the freely available open-source algebra program SINGULAR [2091 is employed. It 
is found that large parts of the parameter space of the NMSSM Higgs potential can be excluded by requiring the 
global minimum to have the electroweak symmetry breaking observed in Nature. This is illustrated by determining 
the allowed and forbidden ranges for some generic parameters of the model. 



6.1. Stationary points in the NMSSM 

The NMSSM Higgs potential is given in Sect. |3.l] with parameters 

''^^K, m|j^ ,rn|^^,m|, Aa, (6.1) 



The quartic terms of the potential (3.8 1 are positive for any non-trivial field configuration, if both A and n are non- 
vanishing. The potential is therefore bounded from below for all cases considered here, and stability need not to be 
checked any further. 

The NMSSM Higgs potential is translated to the formalism described in App. [G] where all Higgs-doublet products 
are replaced by real gauge-invariant functions, Kq, Ki, K2, Ki, and the complex singlet field is decomposed into two 
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real fields according to S = Sr 



In this notation the potential is given by 
1, 



-Re(AK*) {Ki{Sl^ 



o2 



{Kl-Kl-Kl-Kl)^ 

\ml^ {Ko - K,) + ^m%^ {K, + K3) 
+ ms 

-Re(AAA) iK,Sre-K2S,„,) 

+ Im(A^A) {K2Sre + KiS.m) 

+ ^ReMJ(^3^-3^.e5L) 

2 3 2 



(6.2) 



For given values of the potential parameters ( |6.1[ ) we can find all stationary points of the NMSSM by solving the 
systems of equations (Gill, (G13), (G15l) as described in App. [G| 



The initial parameters in the Higgs potential (6.1 1, are translated to the set of parameters as described in Sect. 3.5 



This enables for instance to fix the vacuum expectation values v, tan(/3) and Vg. As a numerical example the parameter 
are chosen as 



A = 0.4, K = 0.3, |^„| = 200 GeV, 
tan/? = 3, rn}j± — 2v, 
signi?«, = -, (5edm = 0, (5^ = 



(6.3) 



and the parameter Vs is varied. The roots of the univariate polynomials are found numerically, with a precision of 100 
digits chosen. The approach allows to use arbitrary precision. The errors of the approximate statements described in 
the following are checked to be under control. For generic values of the parameters 52 complex solutions are found: 
7 corresponding to the unbroken, 38 to the partially broken, and 7 to the fully broken cases. The number of real and 
therefore relevant solutions depends on the specific values of the parameters. 

As expected from the symmetry of the potential, either 1 or 3 solutions sharing the same value of the potential 
are found within the accuracy of the numerical roots. From the computed stationary points only those may be 
accepted as global minima which correspond to the initial vacuum expectation values (up to the complex phases), 
that is which fulfill 



K0 + K3 



tan/3, j2{S?., + SfJ 



(6.4) 



Since for non-vanishing A, k the potential is bounded from below, the stationary point with the lowest value of the 
potential is the global minimum. 

In Fig. |20] the values of the potential at all stationary points for the parameter values (6.3) with varying Vs is 
shown [63 . Each curve in the Figures represents 1- or 3-fold degenerate stationary potential values, where the gauge 
symmetry breaking behavior of the solutions is denoted by different line styles. Excluded parameter regions, where 
the global minimum does not exhibit the required expectation values (6.4 1 are drawn shaded. As is illustrated in the 
figure, substantial regions of the NMSSM parameter space are excluded. For some excluded parameter regions, the 
partially breaking solutions with the required vacuum expectation values (6.4 1 are saddle points. This means they 



can be discarded as global minima without calculation of the other stationary points. However, this is not always the 
case. Obviously from Fig. |20| an upper bound for Vg is found. For the plotted Vg larger than this upper exclusion 



bound the solutions fulfilling (6.4) are still pronounced minima, i.e. the mass matrices have positive eigenvalues, but 
they are no longer the global minima. Also there are non-breaking saddle points with potential values slightly above 
those of the wanted global minimum. It is found that this effect is not coincidental for the initial parameters (6.3 1 
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FIG. 20: Values of the NMSSM potential at its stationary points in dependence on Vs, as presented in W3f . The parameters 
of the Higgs potential are chosen as X = 0.4, re = 0.3, |^k| = 200 GeV, taxi P — 3,mff± — 2?;,signiiK — —,5edm = 5'i^=Q. Each 
line corresponds to 1 or 3 stationary points sharing the same value of the potential. The different line styles denote saddle 
points, maxima, and minima. The labels 'none', 'full', and 'partial' denote solutions of the classes with unbroken (Gil), fully 
broken (G13), and partially broken , G15) SU{2)l x U{1)y, respectively. For solutions of the partially broken class, it is also 
denoted whether they correspond to the 'required VEVs' Vu,Vd,Vs or to 'other VEVs'. Excluded parameter values, where the 
global minimum does not exhibit the required vacuum expectation values, are drawn shaded. 



chosen, but rather a generic feature of the NMSSM. Within the CP conserving parameter range 



Ae]0,l], Ke]0, 1], A„ e ±]0, 2500] GeV, 
tan P € ]0, 50], G ]0, 5000] GeV, 
mH± e] 0,2500] GeV 



(6.5) 



samples are selected, producing the wanted global minimum and typically non-breaking saddle points are found, where 
the relative separation of the potential values for the saddle points and the global minimum is below the per-mill level, 
in many cases even far below. No fully breaking global minima are found for scenarios in the range (6.5 1 where the 
solutions with the required vacuum expectation values (6.4 1 are local minima. Eventually, there are examples, where 
CP conserving parameters with the "wrong" global minimum produce the wanted global minimum if a non-vanishing 
phase S'^ is introduced. 
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7. SUMMARY 

The MSSM suffers from the /z-problem, that is the dimensionful /ti-parameter has to be adjusted by hand to the 

clcctrowcak scale. In the NMSSM an efFcctivc /i-tcrni is generated dynamically. To this purpose an extra singlet 
superfield is introduced. Via spontaneous symmetry brealcing the Higgs-boson singlet acquires a vacuum-expectation- 
value and generates the required /U-term dynamically. 

However, as we have seen, this additional term imposes a continuous symmetry, the Peccei-Quinn-symmetry 
in the superpotential. An additional cubic self-coupling term of the singlet is introduced in order to break the 
Peccei-Quinn-symmetry which would lead to an unobserved axion. That is, the PQ-symmetry is promoted to a 
Ks-symmetry, which, since spontaneously broken, leads to the formation of dangerous domain walls. The loophole is 
to impose higher order operators and additional symmetries, breaking the Ks-symmetry but not disturbing physics 
at the electroweak scale. 

In the NMSSM we encounter two more Higgs bosons as well as a fifth neutralino; compared to the MSSM. The 
modified mixing matrices were recalled in the Higgs sector as well as in the neutralino sector. Also the physical 
potential was derived, yielding, in particular, additional trilinear ^-parameter terms. 

The theoretical and experimental constraints on the model were discussed. On the theoretical side restrictions arise 
from the symmetries of the physical potential. Thus a global minimum which is electric- or color-charge breaking is 
forbidden, as well as a minimum with the wrong electroweak symmetry breaking behavior. Stability of the potential is 
of no concern, since for a non-vanishing Peccei-Quinn-symmetry breaking K-parameter quartic terms in the potential 
ensure a potential which is bounded from below. 

Further theoretical constraints come from the requirement of perturbativity of the couplings up to the GUT or 
Planck scale. Of course, this condition relies on the viability of the model up to large scales. In a series of recent 
publications the argument of fine-tuning got a lot of attention. On a quantitative basis, studies of fine-tuning in 
different scenarios in the NMSSM as well as compared to the MSSM were performed, revealing that in general there 
is much less fine-tuning in the NMSSM. This fine-tuning mostly arises from the stop-loops entering logarithmically 
in the Higgs-boson selfenergy. In particular, the LEP bounds require large radiative corrections accompanied with 
very large stop masses in the MSSM. Compared to the MSSM, the Higgs sector is much less constraint, that is, the 
NMSSM passes the LEP constraints with much less fine-tuning. 

We proceeded with a discussion of the experimental constraints. First of all the electroweak precision measure- 
ments were discussed, which agree well with the SM predictions. Especially the invisible Z-boson decay puts severe 
constraints on the NMSSM, but also the constraint of limits on neutralino/chargino- and charged Higgs-boson pair 
production at LEP was discussed. 

The impact of the muon anomalous magnetic moment, measured at BNL, on the NMSSM was discussed. From 
chargino-sneutrino and neutralino-smuon loops we get additional contributions to the anomalous magnetic moment 
compared to the SM. Based on clear SM predictions which axe not yet available strong constraints on the parameter 

space in the NMSSM could be derived. 

Also the 6 ^ S7 decay was considered, which, since loop-induced, is very sensitive to new particles, which couple 
to SM particles. Here the measurement agrees well with the SM prediction, yielding an additional constraint for new 
contributions which arise in extensions of the SM. Also in this respect the NMSSM does not violate this constraint 
significantly in large parts of parameter space. 

The new five-year WMAP data give accurate predictions of cold dark matter in the Universe based on the ACDM. 
Since supersymmetric models, respecting matter-parity, predict a new stable particle, they provide a natural cold 
dark matter candidate, the LSP. The WMAP constraint on the LSP turns out to be very strong in the NMSSM and 
restricts the available parameter space enormously. In the constraint NMSSM, where unification of the scalar mass, 
of the gaugino mass, and of the trilinear A-parameter at the GUT scale is assumed, the theoretical and experimental 
constraints are shown to be highly restrictive. 

Strong first order electroweak phase transitions in order to account for the observed baryon-antibaryon asymmetry 
in our Universe require a cubic term in the physical potential. With this respect in the SM as well as in the MSSM 
this mechanism relies on generically small loop contributions in the effective potential. In contrast, in the NMSSM, 
the additional trilinear terms allow to accomplish for this mechanism of baryogenesis without large fine-tuning. 
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Some recent parameter scans were reviewed with quite interesting results, constraining the NMSSM parameter 
space. Some emphasis was placed on a "no-lose" -theorem, that is, the question whether at least one Higgs-boson 
is detectable at the LHC with high integrated luminosity. The Higgs-to-Higgs decays are found to be rather 
difficult to detect, which however correspond to large parameter space. New ideas to detect such signatures were 
reviewed. However there seems to be some parameter space left, where all supersymmetric partner particles could 
escape detection. An electron-positron collider could close this gap in detecting the recoil mass in signatures of 
Higgs-boson production in invisible Z-strahlung. Some publications were reviewed discussing how to distinguish 
the NMSSM from the MSSM. The fifth neutralino may lead to a very different signature in colliders. In case the 
LSP is a singlino-like neutralino, this LSP has only suppressed couplings to non-Higgs particles. Since eventu- 
ally all superpartner particles decay into the LSP this would cause very different signatures, possibly displaced vertices. 

Eventually, we draw attention to the determination of the global minimum in the Higgs potential. An algebraic 
method bases on Groebner bases computation was introduced, showing a surprising rich structure of stationary 
solutions. However, this approach is yet limited to studies of the tree-level potential. 

The NMSSM is an intriguing model, which deserves a lot of attention, since it is a c;ohcrent supersymmetric 
extension of the Standard Model. Moreover the NMSSM complies with collider and cosmological precision data and 
could be discovered at the LHC. 
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APPENDIX A: CONVENTIONS AND ABBREVIATIONS 



We use the space-time metric (g^jy) = cliag(l, —1, —1, —1). 
The generaUzed Pauli matrices are 

'^0 = CTO = [1 l) , '^l^-'^l^ (? jj' '^2 = -^2 = r- ~q) , ^3 = -^3=(J_M- (Al) 



The e symbol matrix is 



We use the usual representation for the Dirac matrices 



1 
-1 



Further the projectors are 





and we use the definition 



(A2) 



7.= (,yo> and 75^(;^_«g. (A3) 



= " V j ' ^« = - (0 I2 ■ 



The conventions of the vacuum expectation include the square-root of two: 

{H,} = H/^) , = e^*" f ° ) , (S) = e''^'v,/V2, (A5) 



ta.n{P) = tf3=—, with v= Jvl + vl^i 2^6 GeY . (A6) 

Wc use also cp = cos(/3), sp = sin(/3), cot^ = cot(/3). After electroweak symmetry breaking we thus get the tree 
level gauge-boson masses niw = g2v/2 and mz = mw/cw, where we use the abbreviations for the Weinberg angle 
cw = cos(^vr) as well as sw = sin(^vr)- 
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APPENDIX B: COMPUTER TOOLS FOR THE NMSSM 

Here we list some frequently used computer tools for the calculations in the NMSSM. In the parameter scans 
discussed in Sec. |5.4| these tools appear in different contexts. 

• NMHDECAY [1071 1108] This Fortran code computes the masses of all sparticles and Higgs bosons. The couplings 
and decay width of the Higgs bosons are also calculated. The input parameters are A, k, A\, A^, tan(/3), 
/i = Xvs taken at the electroweak scale. The computation of the Higgs-boson spectrum is done including leading 
electroweak corrections as well as certain two loop terms. The decay width refers to HDECAY |210j . but without 
taking into account three body decays. 

Current experimental exclusion limits are taken into account. The program code may be downloaded from the 
url: http : //www . th . u-psud . f r/NMHDECAY/nms smtools . html| 

• NMSPEC PTT] In contrast to NMHDECAY here the soft-breaking parameters have to be specified at the GUT 
scale. Also the spectrum and couplings are computed from this input. Download at url: http : //www . thT] 



u-psud ■ f r/mHDECAY/nmssmtools . html 



• MicrOMEGAS |212l I213j This Tool calculates the relic density of a stable massive particle together with the 
rates for direct and indirect detection of dark matter. MicrOMEGAS includes already a model file for the 
NMSSM in addition to various other models. It may even be extended to further models by the user. Download 



available at url: |http : //wwwla pp . in2p3 . f r/lapth/micromegas 



• CompHEP [214' This package calculates total and differential cross sections at tree-level accuracy. Multi-particle 
final states in collisions as well as decay processes can be computed in a completely automatic way. We refer to 



the url: http://coiiiphep.sinp.msu.ru 



LanHEP |215j LanHEP computes the Feynman rules for a given Lagrangian. The initial Lagrangian can be 
written in a compact form. The output Feynman rules are given in terms of physical fields and independent 
parameters. It can also be used to directly generate a model file as input for CompHEP. Download at the url: 



http : //theory . sinp . msu . ru/'^semenov/lanhep . html 



APPENDIX C: CONSTRUCTION OF A SUPERSYMMETRIC MODEL 

Before we start to present the Lagrangian of the NMSSM we want to sketch how an arbitrary supersymmetric 
model is constructed. Based on this sketch it is easy to construct the specific NMSSM Lagrangian in the next section. 
Here we follow closely the excellent introduction given in Ref . [21] . First we recall the meaning of matter parity. 

1. Matter parity 



Writing down the superpotential in the MSSM (2.2 1 or the NMSSM ( |3.2| , lepton and baryon number violating terms 
are omitted. Additional lepton number (B) and baryon number (L) terms in the superpotential with dimensionless 
couplings would be 

Wab, al =- ^Ae (L'^eL) + X'd{L'^eQ) + ^\"ddu . (CI) 

The non-appearance of such terms in the superpotential can be gained by imposing an additional symmetry principle, 
called matter parity |12l 12161 12171 1218] or R-parity |10j . Matter parity is defined as a multiplicative quantum number 

Pm - (-1)3(^-^) (C2) 

such that only terms in the Lagrangian or the superpotential are allowed with multiplicative matter parity P/\/ = +1. 
In this way the forbidden terms are excluded from the model. Equivalently one can impose also R-parity instead of 
matter parity, defined as the multiplicative quantum number 

with s the spin of the particle. Since the product of (—1)^* is always 1 in angular momentum conserving interaction 
vertices, R-parity is indeed equivalent to matter parity but has the advantage to give for the SM particles and Higgs 
bosons Pr — +1 and for the superpartners Pr = — 1. R-parity conserving immediately translates to the fact that the 
lightest supersymmetric particle (LSP) is stable. Since the NMSSM like the MSSM is R-parity conserving we expect 
to have a LSP which might be a candidate for the up to now missing cold dark matter. 
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2. Chiral supermultiplets 

The boson and fcrmion fields are cast into chiral supermultiplets. Each of the n boson fields of the theory, (jji 
with i — l,...,n is accompanied by a Weyl-fermion ipi and an additional auxiliary field Fi necessary to close the 
supersymmetry algebra off-shell. Note that the auxiliary fields Fi are no physical fields since they do not propagate. 
The supersymmetric Lagrangian, that is the Lagrangian which is invariant under supersymmetry transformations 
which turn bosons into fermions and vice versa, is 

/:chiral. free = - (9^^^^ " + F'* F, . (C4) 

The generalized Pauli matrices are do = = I2, ~ ^fi, (T2 = — cr2, o'l = ^cs- The most general set of 
renormalizable non-gauge interactions of these chiral supermultiplets are 

>Cehiral, int = -^W'^ ^,^P, + W F, + C.C. , (C5) 

where as well as W^^ are determined from one function, the so-called superpotential W: 

W = ^M^'M, + Iv'^'Mjcl^k (C6) 



with 



SW 1 I. 



(C7) 



5^W -u 



This general form of the interactions ( C5 ) is dictated by supersymmetry itself. From the Lagrangians ( C4 1 and ( C5 1 
we find the equation of motion 

F, = -W* , (C8) 
thus the auxiliary fields F^ can be expressed in terms of the scalar fields. 



3. Gauge supermultiplets 

The gauge boson fields A^^ are paired with Weyl fermions A" into gauge supermultiplets, where also auxiliary fields 
are needed in order to close the supersymmetry algebra off-shell. The adjoint representation of the gauge group 
is denoted by the index a here. The Lagrangian of the supersymmetric gauge supermultiplet fields is 

C„ = - - iX^d^^D^X- + ]^D-D'' , (C9) 

where, as usual, the Yang-Mills field strength is — dfj_Af, — d^A^;^ — g f"''^'^ A^^A^ and the covariant derivative of the 
gaugino field reads 

D^X'^ = d,x - gr'^A^x' (CIO) 

with structure constants /"''^ and gauge coupling g. 



4. Interactions 



As usual the gauge interactions of the bosons and fermions are given by turning the partial derivatives in the 
kinetic terms in the Lagrangian (C4| into covariant derivatives. Suppose the chiral supermultiplet fields transform 
under a gauge group in a representation with hermitian matrices satisfying [T", T''] = if'^^T". Then the covariant 
derivatives of the multiplet fields are 



D^(j,,^d^(i3, + igA1{T''4>), 



(Cll) 
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Thus, the gauge invariant form of the Lagrangian ( |C4[ ) becomes 

-Cchirai = -{D^^')^D'^^,) ~ tij^a^^D^i;, + F"F, . (C12) 

Further we have to introduce all additional terms in the Lagrangian, which do not violate any symmetries of the 
theory. The following Yukawa terms are neither forbidden by gauge invariance nor by renormalizability. The couplings 
of these Yukawa interaction terms are determined by supersymmetry. The Yukawa couplings are 

/:Yukawa = ~ ^5 [{<l>*T^i^')\^ + X'^H^^ T'^ <l>)] + 9{<jy*T'^ c^)D- ■ (C13) 

From £gaugc and /^Yukawa we get the equation of motion for the auxiliary field 

= -g(0*T"0) (C14) 

and we see that like for the auxiliary fields Fi we can express Z?" in terms of scalar fields. We emphasize that the 
Yukawa coupling of a fermion with a scalar and a gaugino is determined by the gauge coupling, one of the firm 
predictions of supersymmetry. 



5. Soft breaking terms 

We know from experiment that supersymmetry must be broken. This breaking is expected to be spontaneous, 
such that the Lagrangian is invariant under supersymmetry transformations but the vacuum is not. The breaking 
mechanism itself is up to date unknown. In order to keep the theory as general as possible all explicit breaking 
terms are introduced which do not lead to quadratical divergences. A necessary condition for this is to introduce only 
Lagrangian terms with couplings of positive mass dimension. The most general soft breaking terms are [219^ 

Csoft = -(m2);.0J - i (MaA'^A'' + c.c.) - Qfo'V'^V,,^^- + ^a*^''=A*^'V.'/',0fe + c.c}j . (C15) 

The first term gives masses to the scalar superpartners and the second one masses to the gauginos. In this way the 
degeneracy among the superpartners is removed. The bilinear (b^^) and trilinear (a*-''^) terms in the bracket are terms 
associated to the superpotential bilinear and trilinear terms. 



6. Complete supersymmetric Lagrangian 



Eventually we have the general supersymmetric Lagrangian of a renormalizable theory consisting of the chiral- and 
gauge- supermultiplets in gauge invariant form as well as the Yukawa and the soft breaking part: 



^ -^chiral ^" -^gaugc ^" -^Yukawa "t" -^soft ■ 



(C16) 



We may also isolate the potential part of this Lagrangian. We have the F term of £chirai, the D term of >Cgaugc and 
the scalar part of >Csoft, altogether 



V = Vf + Vd + Vsott 
= ~F"F, 



(C17) 



where the equation of motion for the auxiliary fields (C8l and (C14l have to be inserted. 
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APPENDIX D: FEYNMAN RULES OF THE NMSSM 



In this appendix some of the essential Feynman rules of the NMSSM are derived. With the previous work in 
appendix |C] this is straightforward. We focus on the Higgs sector and restrict ourselves to the trilinear couplings. 
In Ref. [16 the Feynman rules for a generic supersymmetric two-Higgs-doublet model are derived, taking also an 
additional singlet into account. But this work is restricted to the case that the singlet does not mix with the other 
neutral Higgs bosons - in contrast to the case of the NMSSM. The Feynman rules for the NMSSM can also be found 
in Refs. [571 1107] . Here we want to present the Feynman rules for the general case, including CP violating complex 
phases in the NMSSM parameters. In this case, the scalar Higgs bosons mix with the pseudoscalar ones, as described 



in section 3.3 Note that we perform the P rotation of the Higgs-boson doublets explicitely. 

Arrows on the field lines in the Feynman diagrams denote the charge flow. In case the Feynman rules depend on 
a momentum of a field, the direction of momentum is denoted by an extra arrow. Further, we use the abbreviations 
for the Weinberg angles, sw = sm{9w), cw = cos(6'n'), /3 is the mixing angle of the Higgs doublets with tan(/3) = 



The mixing of the Higgs bosons is determined by the mixing matrix (Rij) and 



tfj = Vu/vd and w = < + 

originates from the diagonalization (3.22) with indices i,j = 1,...,5. Since the resulting entries in this matrix are 
very involved it is advantageous to perform the diagonalization in a numerical way. We use the short notation 

{F G) = F{d^G) - {d^F)G. 



1. Higgs-boson gauge-boson interaction 



The interaction of the Hi ggs b osons with the gauge bosons is determined completely by the gauge invariant kinetic 
terms of the Higgs bosons (C12l 



/:Higg«, kin = ^{D^Hd)\D^Hd) - {D^Hu)HD^Hu) ~ {d^S*){d''S) 



(Dl) 



Note that the Higgs singlet S is supposed not to have any gauge interactions, that is the covariant derivative in this 
case is just the usual partial derivative. The covariant derivative for a field with electroweak interactions reads (Cll I 



igiB^{x)Yw + ig2W;{x) 



(D2) 



with hypercharge operator Y^k and isoweak operators Iw and the corresponding gauge fields B^{x) respectively 
W^{x) . The Higgs doublets have weak isospin and hypercharges Iw Hd — cf"" /'^■Hd and Iw Hu = cr°/2i?„ 
with a = 1,2,3 and Yy/Hd — —l/2Hd and Y^Hu = +l/2-ffu (see Tab. [T|). T he Higgs bosons are singlets 
with respect to the SU{'i)c gauge group. Inserting the covariant derivative (D2| into (Dll and using the pa- 



rameterization (3.9 1 for the Higgs-boson doublets in the unitary gauge we find the following Lagrangian terms 



and corresponding Feynman rules. The 5 x 5-matrix matrix R is defined in (3.22 1 and is to be determined numerically. 



wwww 



'2\\ 



\r\j\r\r\r\j\/\A 




^{P2 - Pi) [ Ri2Cf3 - RilSp + iRv 
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Hi 





w+ 



w- 



Pi 



wwww 



V2 



2. Higgs-boson fermion interaction 



From the superpotential (3.2) we get via (C5l and (C7l the Higgs-boson fermion interactions (note that there is no 

(D3) 



mass term proportional /i*^ in the NMSSM superpotential) 
Inserting the doublets from Tab. |l]this reads 



(D4) 



Now we translate the Weyl two-component spinors into Dirac four-component spinors. For arbitrary Dirac spinors 
decomposed into Weyl components / = (fLjuV and g = {gL,gR)J we get Jg = f-fog = flgn + f^gL and 
fl59 = -flgn + /jf.9L, that is we have /j^gi = JPl g and flgn = JPr g, where Pl/r = (I4 ± 75)/2. Using the 
parameterization (3.9 1 and four-component Dirac spinors u = {ul^urP' and d = {d^^djip- we get 



Hff 



Vu 

\/2 



Vd 



u{vu + K + iaulbju - ;y|^("<i + hd + iaalbjd + {yuuPLdH^ + yauPRdHJ + c.c.) . (D5) 



A complex phase in the parametrization of 7J„ can be absorbed into the parameter ?/„. The bilinear mass terms of 
this Lagrangian can be used to relate the parameters yd, y« to the fermion masses. We find 



Vu 



•\/2mvi/S/3 ^ ^/2m\YCj3 

The trilinear terms give the neutral Higgs-fermion-fermion interactions with Lagrangian and Feynman rule 

^HOqq = Hi U [ Ri2+ IClsRiA-y^ U - 



gimd 



Hi d\Ri\ + iSf3Ri4'y5 d 



as well as the charged Higgs-fermion-fermion interactions (where we applied the (3 rotation (3.181 and (3.25 1) 
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m u 



(^ucotp + mdtp + {rriucotfs - mdtfj)^^^ d + c.i 



The Lagrangian for the Higgs-lepton interaction may be derived in a quite analogous way. 



(D6) 



(D7) 



(D8) 
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Hi 




92 



Hi 



Ril + is/}R^ij5 



H, 



R 



il 



H+ 



rriucotp + mdtp + {niuCotp - mdtp)"/^ 



3. The NMSSM potential 



In Sect. [6] the potential of a general supersymmetric theory was given. Here we want to derive the full potential 
explicitly for the NMSSM. We get contributions to the potential form the F-terms, £)-terms, as well as from the soft 
supersymmetry breaking terms £soft • The i^-terms read with ( C8 1 



Vf ^ F*F' ^ W*W' = 



SW 



(D9) 



where are all scalar fields occurring in the chiral supermultiplets, that is with view on Tab. |l] = {ul, d-L, 
d*j^, Ve, CL, e|j, H^, H^, H+, S). In contrast to the MSSM the NMSSM has no dimcnsionful coupling in the 
superpotential. As already mentioned this means that the function Wi in ( C7 1 has no term proportional to /i'-' . After 
writing the isospin products out, the superpotential becomes 

W = u*^yu{uLHl - dLH+) - dMuLH- - dLH^i) - ~e*„ye{^eH- - giiJ^) + \S{H-H+ - H^dK) + \^S^ (DIO) 
and we can construct the functions Wi and arrive at 



Vf 



\X{HueHd) + kS^ 



UuiQeHu) + VdiQeHd) + ydLeHd) 



ydd*RH^d - yuU*RH+ + lyee^H^l + \yee*RH^ [ ■ 



\yuURUL - \SH^ 



(Dll) 



ydd*RdL + yee*ReL- XSH" + \SH+ - yddRUL - yee*jii>, 



XSH^ ~ yuU*RdL 
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The next contribution to the potential conies from the D-terms in (C17l. With the equation of motion for the 



auxiliary field Z?" in (C14l we get 



(D12) 



where 4)i denotes again all scalars of the chiral supermultiplets and ga is the gauge coupling corresponding to the 
gauge group generators T" and there is an implicit sum over the adjoint index a. For the SU(2)l gauge group the 
gauge coupling is 52 and the generators are half the Pauli matrices; T° = cr"/2 with a = 1,2, 3. With view on Tab. [l] 
we have the scalar weak isodoublets cf) — {Q, L, Hu, Hd) available and get for instance for the doublets Q and L a 
contribution to (D12l 



(D13) 



where r, s,t,u — 1,2 and we used the identity <j!^s(Tf^ — 25ru5st — 5rs5tu- For the U{1) gauge group the gauge couphng 
is gi and the generator is the hypercharge of the fields T" — Yw given in Tab. |lj We have the following scalar fields 
with non-vanishing hypercharge (p = (Q, w'^, d*j^, L, e^, -f/„, Hj). Exemplarily, for the Higgs doublets Q and L with 
hypercharges Uq and y^, respectively, we get a contribution to (D12) 



Igl yi (Q^QKDl) 



(D14) 



Collecting all D-term contributions we arrive at 



Vd 



9l 



4|ijtQ|2 „ 2(ijti7„)(gtQ) + 4|^t£|2 _ 2{HIH„){Dl) + 4|gtL|2 - 2{Q^Q){Dl)] 



(D15) 



2 V6 



>UrUr 



P*ReR 



Now we come to the explicit soft supersymmetry breaking terms, that is terms which do not introduce divergences 
at higher order in perturbation theory. Writing down all terms as given in (C15l we have 



•Csoft — 



t , 



XA^SiH^eHd) + -A^S' + c.c. 



'pi - ml\uR\'^ - TO^~MflP - ml\eR\ 



u*RyuAu{QeHu) - d*j^ydAd{QeHd) - e*RyeAe{LeHd) 
MiBB + M2WW + Mggg + c.c. 



c.c. 



(D16) 



with new soft parameters rnj^^, rrig, mg, 



f^L, ml, m^, m% 



Ax, A^, Au, Ad, Ae, Ml, M2, Mg. The idea is. 



that these parameters one day emerge from the supersymmetry breaking mechanism. As long as this mechanism is 
unknown, this supersymmetry breaking soft terms parameterize the generic model. Note that due to missing bilinear 
terms in the NMSSM superpotential there are no associated soft terms of this kind in the corresponding soft terms. 

The potential is given by the scalar part of £so/t, that is without the last line on the right hand side and with a 
negative sign since it is a potential. The Higgs-boson part of the soft breaking terms is given by the first line on the 
right hand side of (DI6I. 



55 



4. Higgs sfermion interaction 



The Higgs sfermion interactions follow from the F-terms (Dill, the D-terms (D15I as well as from the soft 
terms (D16). The corresponding Lagrangian is 



i^Hqq — — 



VuiQeHu) - VdiQeHd) - yuU*iiHl - ydd^H^ - ydd*RH^ - yuU*i^H+ 
yu{uBulXSH^ + URdlXSH- + c.c) + yd{dRdlXSH° + dRulXSH+ + c.c) 



9l_ 
12 



{HlH^){Q^Q)~{HlHd){Q^Q) 



(D17) 



9l 



u*RyuAu{QeHu) - d*RydAd{QeHd) - e\yeAe{LeHj) + c.c 



The Higgs-slepton interaction is derived in a quite analogous way and not given here explicitly. The Feynman rules 
follow from the parameterization (3.9 1, the rotation of the Higgs gauge eigenstates into mass eigenstates (3.221 and 



the mixing of the and qji gauge eigenstates into mass eigenstates with mixing matrix Rq, 



cos{0) sin(0)\ (qA _ 



-sin{e) cos{e) \qii 



Rfj 



(D18) 



that is by the diagonalization of the mass matrix M; diag(mgj,TOq, ) = RqMR^ 



5. Higgs-boson neutralino/chargino interaction 

We start with collecting the bilinear neutralino/chargino interaction terms as well as the Higgs- neutralino/chargino 
interactions. 



From the soft breaking terms (D16 1 we get mass terms of the bino and the bino with mass parameters Mi and M2 
respectively. 



^soit,B/w = -\miBB - ^M^W'^W^ + c.c. = -hd^BB - M2W-W+ - ^M^W^W^ + c.c. 



(D19) 



where we used the definition 



V2 



in an analogous way to the vector bosons. 

The chiral interactions ( C5 ) give us additional Higgs-Higgsino terms. 



-chiraLHff 



Insertion of the Higgs-boson doublets and Higgsino doublets from Tab. |T] gives 
1 



(D20) 



1 1 

'2WW-^''^' + c.c. = - - ( X{H'^eHd)S + X{H^eHd)S + X{H^eHd)S + nSSS ) + c.c. (D21) 



^chiraLHH = " 3 ( S - XH^^H^J XH",H^J + XH' H+ S A5i/0ff„" + XSHJH+ + ^SSS ] + c.c. (D22) 



Further Higgs-Higgsino interactions come from the Yukawa interaction ( |C13 1 



-Yukawa, 



= - \/2.g(0|T'^V')A'^ + c.c. = 

- V2g2 (nl^Hd + Hi^H^ W'^ - V2gJ - ]-HlHd + IhIH^ ] Bo + c.c. 



(D23) 
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We get, using p20 l 
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(D24) 



Now, since we have collected all necessary Lagrangians we can easily give the neutralino mass matrix, as presented 
in (4.1 1. We choose the basis '0'' — {B, W^, H^, H^, S). In this basis the neutralino mass part of the Lagrangian reads 



M^oiIj" + C.C. 



(D25) 



with 



Ml —C/iswfnz spSwiTLz 

Ah cpcwmz -Sf3Cwmz 

-cpSwiTT-z CfjCwmz —Xvg/V^ — Av„/\/2 

s/jswrnz —spcw'mz —\vs/V^ —\vd/V^ 

-Xvu/V2 -Xvd/V2 



(D26) 



The upper left 2x2 block arises from the soft breaking mass parameters of the bino B respectively wino W (D19l. 
The lower right 3 x 3 blo ck, mixing the Higgsinos among themselves, comes from (D22l and the remaining part from 
the Yukawa terms (D24l, where we get bilinear terms from the vacuum conditions {H^) — Vd/V^, (ff") = Vu/V^ and 
(5) = Vs/V^. The mass eigenstates follow from a unitary rotation, that is 



Xi = Ujiip^ with diag(m^o,m^o,m^o,m^o,m^o) = U*M^oU^ 



(D27) 



with i,j — 1, 5. Conventionally the neutralino masses are defined to be in ascending order, that is Xi is the lightest 
neutralino. The result of this diagonalization is rather involved and in practice performed numerically. 

From the chiral interactions (D22l and the Yukawa interactions (D24) we get the neutral Higgs boson neutralino 
interaction 



- AtLdtl.^b - AtL^tldb - AbHdtL^ 



kSSS 



^[{H'iYH'dW'-{KYHlW^ 
^^[{H'.YHlB'-iKrHlB' 



(D28) 



Using the Higgs boson parameterization ( 3.9 ) , diagonalizing the Higgs bosons ( 3. 18 1 , ( 3.22 1 , and the neutralinos ( D27 1 
we have 



Ch<,«x-=A'"'H,x'IxI + c.c. 



with 



- X{{Rn + iR,ASp)U*M5 - {R^2 + ^R^iCp)U*:,Ul, + (i?,3 + lR^^){U*^Ul^ + U*,Ul,))] 
Now we write the products of Weyl-fermions in terms of Dirac four component spinors 



with i = 1 , . . . , 5 . 



(D29) 



(D30) 



(D31) 



that is we have for the products in terms of Dirac four-spinors 



X°xl = X^PlxI and (x°X°)* = X°7oPl7oX? = xIPrX° = x"PrxI 



(D32) 
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where the fields denoted with a tilde symbol are the Dirac spinors and the last equation holds only for Majorana 
spinors [16] . With ( D32 1 we get 



(D33) 



From the soft terms Lagrangian (D19), and the chiral terms Lagrangian (D22l we can immediately construct the 
chargino mass matrix. With the chosen basis ij}'^ = (PF"'", iJ^, iJ^) the chargino mass part of the Lagrangian 
reads 

1 

2< 



-Cchargino mass = -^{^^)'^ M^±^^ + C.C. 



(D34) 



with 



The unitary rotations 



xt 



can be used to get the mass eigenstates 





X 



X 



M2 ^/2s|3mw\ 
Afs/V2 ) 



W+ 
H+ 



u- 



w- 



diag(m^±,m^±) = (;7+)*X(C/-)t 



Explicit diagonalization yields the chargino masses 



m ±, ± 

Xi /Xa 



^ (^|Af2p + |Ap«2 + 2ml. T v/(|A/2|2 + |A|2«2 + 2m2^)2 - ilXv^Ah - sin(2/3)P 



(D35) 

(D36) 
(D37) 
(D38) 



From the chiral interactions (D22) and the Yukawa interactions (D24) we get the neutral Higgs boson interaction 
with two charginos. 



-Cffx+x- = \SH~H+ - g2{H° *H~W+ + Hi *H+W-) + c.c. 



(D39) 



Using the Higgs-boson parameterization (3.9), the rotation of the gauge eigenstates into mass eigenstates of the Higgs 
bosons (3.181, (3.221 and of the charginos (D36l, we get the trilinear terms 



C 



,Xj xt + c.c. 



(D40) 



with 



Transferring to four-component Dirac spinor notation: 



xt ■■- 



(D42) 



where the charginos denoted with a tilde symbol are the Dirac four-spinors and the charginos without a tilde symbol 
the two-component Weyl spinors. This means we have the translation xt ~ ^Lxt ^'^^ X4 — PRxt ' that is x7 — 
xt Pl and arrive at 



Xi 



Chx+x^ = B'^^Hd^PLxt + c-c. = xt {B^'^Pl + {B^'^Y Pr) xtH^ 



(D43) 



Finally we have the charged Higgs interaction with one neutralino and one chargino. Collecting the appropriate 
terms in (|D22| and (|D24| we find 



A 



(D44) 



.91 



With the parameterization 



nos (D27l and the charginos 



58 

3.9 1, the P rotation of the charged Higgs bosons (3.251, the rotation of the neutrah- 

(D45) 



D27 1 we have 



where we defined 



A 

2V2 



{spiu+)i^u;,-c0{U')k2Uj5) 



+ ^{cp{U+)l^U*^~sp{U-)u2U,2) 
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In Dirac four-component notation we find eventually 



(D46) 



(D47) 




6. Higgs-boson self interaction 



The Higgs boson self interactions follow from the potential, that is i^-terms (Dill, the _D-terms (D15l as well as 
from the soft terms (DI6I. Here we present the trilinear part of the corresponding Feynman rules |107j . where it is 
defined 



Pijk — PilPjmRkm ^" RilRkniPjm 



Pijk — Ril{RjmRkm ^" RkmRj^n) j 



(D48) 



with R given in (3.221. 



^ p233 I p244 2 I p255\ 

. pl33 

ijk ~ ijk 

.- . p311 , p344 2 , p344 2n 



ijfe 

pl44 2 , pl55^. 



Xk 



/ p313 pl55 
'^u\^ijk ^ijk 



>323 p255 , 9p345„^\ 
ijk ~ -'ijk 'T^^ijkC-f}) 

O,, ^p213 p344„ „ , p245„ , pl45^ ^ 



p355 
ijk 



3%/2 



p333 
^ijk 



p355 



.91 + g| 



p213 



p245„ pl45„ p344^ 



ijk ^ijk ^ ^ijk '^ff ^ijk ''(3) 
— Vd{Pij'k — Pijk ^ Pi^k'^13 + Pifk^p) 



VuiP^ 



p211 



yUsRiS ~ VuRilCpSp - VdRi2C0Sj3 

+ 2KXVsRi3Ci3Sf3 + V2XAxRi3Cf3Si3 
+ -J ( Vu{Ri2c} - Ri2sl) + Vd{RilSp - RiiC^) 

+ ^ ( Vn{Rz2 + '^RilCfjSp) + Vd{Rtl + 2Ri2CpSj3) 
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APPENDIX E: THE CONSTRAINED NMSSM 

In the NMSSM there is a very rich parameter space, mostly originating from the soft-breaking terms. In contrast, the 
remaining part of the model is highly restricted by supersymmetry. Essentially, the soft-breaking part parameterizes 
the missing knowledge of the supersymmetry breaking mechanism. Motivated by the unification of the gauge couplings 
at the GUT scale, a rather generic approach is to assume unifications at this scale of the soft-breaking parameters. 
The NMSSM with this unification assumption is called constrained NMSSM, or cNMSSM, introduced analogously to 
the constraint MSSM, or cMSSM. The initial set of parameters is reduced by imposing soft parameter unification at 
the GUT scale of about Aq = 10^^ GeV. Explicitely, the parameters, taken at the GUT scale are: 
Trilinear >l-parameter unification, 

^0 = Aa(Ag) = A„(Ag) = ^„(Ag) = AdiAc) = A,{Ag) . (El) 

Gaugino mass unification, 

Ml/2 = Mi(Ag) = M2(Ag) = M^(Ag) . (E2) 

Scalar mass unification, 

ml = ml^{AG) - m|,JAG) - m|(AG) = m|(AG) = mi(AG) = mi(AG) - m}{AG) = m?(AG) . (E3) 

With these unification assumptions there remain only five undetermined parameters in the constrained NMSSM which 
may for instance be chosen as the dimensionless parameters [220] : 



VtiAa), A(Ag), k(Ag), ™o/Mi/2, ^/Afi/2 
Some studies of the cNMSSM are presented in Sect. |5.4[ 



(E4) 
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APPENDIX F: FIELD-DEPENDENT MASSES 



In this appendix the field-dependent masses of the top, stop and of the gauge bosons are given, taken from [62 . This 
field-dependent masses arise from the Lagrangian before spontaneous symmetry breaking occurs. All field-dependent 
masses are denoted with a bar sign in order to distinguish them from the physical masses. 

The quark masses are expressed in terms of the Higgs fields and vacuum expectation values, 



(Fl) 
(F2) 



where, recalling the neutral part of (3.9) 



i/o = {vd + hd + iad)lV2 , Hi - e*^" {v,, + K + ia^)/^ , S = e**= (w, + hs + ia,)/\/2 . (F3) 

The physical masses of the particles arise, if we take the vacuum-expectation-values of the Higgs bosons, (H^) = v^/V^, 
{Hi) = e^-vj^, (S) = e''^'vjV2, and we get 



/ - 2\ 2 1 I |2 2 / - 2\ 2 1112 2 

(™6) =rnh^ -\yb\ v^, (mj) = = -\yt\ 



The field dependent masses of the gauge bosons are 

\ttO\2 , |TjO|2 



The gauge bosons masses are 



2 _ + ,2 



/2i 2\ I - 2 \ 2 32/2, 2\ 



(F4) 



(F5) 



(F6) 



The field-dependent top- and bottom-squark masses are 



m: 



_ 1 
■ti,2 ~ 2 



+ < + ^^(I^SP - \Hl\') + 2\y,mi\' 



±J (m? - ml + \ Ll - -gl \ - \Hl\^)Y + 4|y,P|A^i/S - A\Hl* 



(F7) 



mr — - 

bi.2 2 



1 



bR 4 

The masses of the squarks at the vacuum are 



±J {ml -ml--( gl - -gl ) - + A\y,\^\\SHl - AlH^* 



(F8) 



(m? ) 

^ tl,2' 



2 1 

"^i^.^ = 2 



m.~, + mr + 

1 ta 



52 + 5l / 2 



ivi-vl) + \yt\-'vl 



±Jml - ml + ^^-^^{vj - vl) + 2|ytP((l|APt;2 _ R^y^ t&nl3)vj + - R,v, cot (3)vl) 



, (F9) 



\ 61,2' 



fcl,2 2 



2 I 2 

mr, + mr 
1 br. 



52+5l/2 2\ , \ |22 

5 — {Vd-vj + m Vd 



±\ m? -m? - ^^-^{vj - vl) + 2|yb|2((i|A|2z;2 - R,Vs cot (3)vl + - RbVs tan/3)t;2) 



with the abbreviations 



R, = ^Re(A^,e'(^+^)), I, = -^Im(AA,e^(«+^)), (g = &) 



(FIO) 



(Fll) 
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APPENDIX G: GAUGE INVARIANT FUNCTIONS 



In this appendix let us discuss the introduction of gauge invariant functions which turn out to be very helpful! 
in studies of two-Higgs-doublet models and further extensions of the SM. In this we follow closely [63 l 1205] . Let us 
consider the tree-level Higgs potential of models having SU{2)l x U{1)y (weak isospin and hypercharge) electroweak 
gauge symmetry. We draw the attention to models with two Higgs doublets and n additional real Higgs isospin and 
hypercharge singlets. This includes in particular THDMs, where we have no additional Higgs singlets, and the NMSSM 
with one additional complex Higgs singlet corresponding to two real singlets. It is assumed that both doublets carry 
hypercharge y = +1/2 and the complex doublet fields are denoted by 

For the singlets real fields are assumed, denoted by 

^x), i = l,...,n. (G2) 

We remark that in supersymmetric models like the NMSSM the two Higgs doublets Hd, carry hypercharges 
y = — 1/2 and y = +1/2, respectively. This can be translated to the convention used here by setting 



where e is defined in (3.3). Complex singlet fields are embedded in the notation by treating the real and imaginary 
parts of the complex singlets as two real singlet degrees of freedom. 

In the most general SU(2)l x U{1)y gauge invariant Higgs potential with the field content described above, the 
doublet degrees of freedom enter solely via products of the following form: 

with I, j e {1,2}. (G4) 

It is convenient to discuss the properties of the potential such as its stability and its stationary points in terms of 
these gauge invariant quadratic expressions. This was discussed in detail for THDMs and also extended for the case 
of more than two doublets in |205j . We recall the main steps here. 

All possible SU{2)l x U{1)y invariant scalar products are arranged into the hermitian 2x2 matrix 

K := "f^'fA (G5) 

and its decomposition reads 

K,,^\{K^5.,+Kacjt,), (G6) 
where cr° are the Pauli matrices. The four real coefficients in this decomposition are 

Ko = vhr. Ka^{^lvj)Kp a =1,2, 3, (G7) 



where here and in the following summation over repeated indices is understood. The matrix ( G5 ) is positive semi- 
definite, which implies 

Kq > 0, K^- Kl - Kl- Kl>Q. (G8) 



On the other hand, for every hermitian 2x2 matrix (G6), where (G8) holds there exist fields (Gl I satisfying (|G5 



the Higgs doublets. Using the inversion of (G7l 



see |205j . It was also shown in [2051 that the four quantities KQ,Ka satisfying (GBl parameterize the gauge orbits of 



ip\ip,^{Ko + K,)/2, ip\ip2^{K, + iK2)/2, 

+ + (G9j 

4^2 = (Ko - K^)/2, = (Xi - iK2)/2, 

the doublet terms of the potential can be replaced - due to renormalizability - by at most quadratic terms in the real 
functions Kq, Ki, K2, and K3, which simplifies the potential and eliminates all SU{2)l x U{1)y gauge degrees of 
freedom. Eventually, the potential is written in the form V{Ko, Ki, K2, K3, 0i, (/>„). 

To determine the stationary points of the Higgs potential a potential of the form V{Kq, Ki,K2, K3, 0i, . . . , 0„) is 



considered and constraint (G8I is taken into account. The possible cases of stationary points are distinguished with 
respect to the SU{2)l x U{\)y symmetry breaking behavior which a vacuum of this type would have (see ^205j): 
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• unbroken SU{2)l x U{1)y ■ A stationary point with 

Ko = Ki = K2 = Ks^ 0. 



(GIO) 



A global minimum of this type implies vanishing vacuum expectation value for the doublet fields (Gil and 



therefore the trivial behavior with respect to the gauge group. The stationary points of this type are found by 
setting all Higgs-doublet fields (or correspondingly the Kq as well as the Ka fields) in the potential to zero and 
requiring a vanishing gradient with respect to the remaining real fields: 



^)=0. 



(Gil) 



fully broken SU{2)l x U{1)y ■ A stationary point with 

Ko>0, 
Kl - Kl - Kl - 



Ki > 0. 



(G12) 



A global minimum of this type has non- vanishing vacuum expectation values for the charged components of 
the doublets fields in (Gil, thus leads to a fully broken SU{2)l x U{1)y , see |205j . The stationary points of 
this type are found by requiring a vanishing gradient with respect to all singlet fields and all gauge invariant 
functions: 



WV{Ko,Ki,K2,K3,cl)i,. 



(G13) 



The constraints (G12l on the gauge invariant functions must be checked explicitly for the real solutions found. 



partially broken SU{2)l x U{1)y : A stationary point with 

Ko>0, 
K?; - K? - - Kl = Q. 



(G14) 



For a global minimum of this type follows the desired partial breaking of SU{2)l x U{1)y down to U{l)em , 
see [205 . Using the Lagrange method, these stationary points are given by the real solutions of the system of 
equations 



V[F(i^o,^l,if2,i^3,'/'l,.-.,0n) 



-u ■ {Kl - Kl - Kl 



Ki 



0, 



(G15) 



K'. 



Kf - Ki - Ki = 0, 



where m is a Lagrange multiplier, 
for (IgTsI. 



The inequality in (G14l must be checked explicitly for the solutions found 



For a potential which is bounded from below, the global minima will be among these stationary points. Solving 
the systems of equations (Gill, (G13l, and (G15l, and inserting the solutions in the potential, therefore the global 
minima can be identified as those solutions which have the lowest value of the potential. Note that in general there 
can be more than one global minimum point. 



From the mathematical point of view with ( Gl ll, (G13l, (G15) non-linear, multivariate, inhomogeneous systems of 



polynomial equations of third order have to be solved. It is demonstrated that this is possible, even if the number of 
fields is large, like in the NMSSM. The most involved case is given by (G15l, which for the NMSSM consists of seven 



equations in seven indeterminates, namely six real fields and one Lagrange multiplier. In the following an algorithmic 
method to solve (Gill, (G13l, (G15) is described for the case that the number of complex solutions is finite. The 



latter is indeed fulfilled for the NMSSM with generic values for the parameters, and it is automatically checked by 
the method. Note that the gauge invariant functions avoid spurious continuous sets of complex solutions, which are 
found to arise in the case of th e MSSM as well as the NMSSM if the stationarity conditions are formulated with 
respect to the Higgs fields (Gil in a unitary gauge. This is not surprising given the fact, that the gauge invariant 



functions express the contribution of the doublets to the potential by four real degrees of freedom in contrast to the 
five encountered for the doublet components in the unitary gauge. 

The solution of multivariate polynomial systems of equations is the subject of polynomial ideal theory and can 
be obtained algorithmically in the Groebner basis approach |206j . See appendix |h] for a brief introduction to this 
subject. Within this approach the system of equations is transformed into a unique standard form with respect to a 
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specified underlying ordering of the polynomial summands (monomials). This unique standard form of the system of 
equations is given by the corresponding reduced Groebner basis. If the underlying order is the lexicographical order, 
the unique standard form consists of equations with a partial separation in the indeterminates. A variant of the F4 
algorithm |221j is used to compute the Groebner bases. A Groebner basis computation is generally much faster if 
the standard form is computed with respect to total degree orderings and then transformed into a lexicographical 
Groebner basis. The transformation of bases from total degree to lexicographical ordering is done with the help of 
the FGLM algorithm . Finally, the system of equations represented by the lexicographical Groebner basis has to 
be triangularized. The decomposition of the system of equations with a finite number of solutions into triangular sets 
is performed with the algorithm introduced in |223l 1224] . Each triangular system consists of one univariate equation, 
one equation in 2 indeterminates, one equation in 3 indeterminates and so forth. This means that the solutions are 
found by subsequently solving just univariate equations by inserting the solutions from the previous steps. 

The construction of the Groebner basis as well as the triangularization are done algebraically, so no approximations 
are needed. However, the triangular system of equations contains in general polynomials of high order, where the 
zeros cannot be obtained algebraically. Here numerical methods are needed to find the in general complex roots of 
the univariate polynomials. 

In more involved potentials, like the NMSSM, the algorithmic solution is considerably simplified (or even made 
accessible), if the coefficients of the polynomials are given in form of rational numbers. Since rational numbers are 
arbitrarily close to real numbers and moreover the physical parameters are given only with a certain precision this 
does not limit the general applicability of the method in practice. 

All algorithms for the computation of the Groebner basis with respect to a given order of the monomials, the 
change of the underlying order, the triangularization, and the solution of the triangular systems are implemented 
in the SINGULAR program package |209j . The solutions obtained can be easily checked by inserting them into the 
initial system of equations. Moreover, the number of complex solutions, that is the multiplicity of the system, is 
known, so it can easily be checked that no stationary point is missing. 
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APPENDIX H: BUCHBERGER ALGORITHM 



In this appendix we want to sketch the construction of the Buchberger algorithm which transforms a given set of 
polynomials F into a Groebner basis G. Here we follow closely Refs. j63j I208j . The Groebner basis G has exactly 
the same simultaneous zeros as the initial set of polynomials F, but allows better access to the actual calculation of 
these zeros. The general idea is to complete the set F by adjoining differences of polynomials. Before the Buchberger 
algorithm is presented, the two basic ingredients have to be introduced, that is Reduction and the S-polynomial. For 
a more detailed discussion we refer the reader to the literature [20B1 12071 1208j . First of all we recall some definitions. 

Definition. Polynomial Ring 

A Polynomial Ring K[xi, . . . , a;„] = -R'[x] is the set of all n-variate polynomials with variables xi, . . . ,a;„ and coeffi- 
cients in the field K . 

Definition. Generated Ideal 

Let F = {/i, . . . , /„} C K['x\ he finite, F generates an ideal defined by 



I{F) EE ^ E • /. 



In the following we want to consider an explicit example, that is a set F = {fi, f2, fa} C of polynomials 

with rational coefficients: 

/i = 3x'^y + 2xy + y + 9x'^ + 5x - 3, 

/a = 2x^y -xy-y + 6x^ - 2x^ - 3x + 3, (HI) 
/3 = x^y + x^y + ix^ + 2x^. 

The set F generates an ideal /(-F), which is given by the set of sums of /i, /a, and /s, where each polynomial is 
multiplied with another arbitrary polynomial from the ring The summands of the polynomial are denoted as 

monomials and each monomial is the product of a coefficient and a power product. 

Further, we introduce an ordering (>-) of the monomials. In the lexicographical ordering (>-iex) the monomials are 
ordered with respect to the power of each variable subsequently. The ring notation Q,[x, y] defines y )^iox x, that is for 
the lexicographical ordering of monomials powers of y are considered first, then powers of x. Explicitly, this means 
2x^y^ >-\cyi hxy^ because the power of y is larger in the first monomial and 2xy'^ '^Xck ^V^, because both monomials 



have the same power of y, but the first monomial has a larger power of x. The monomials of the polynomials (HI ) from 
the ring (Q[a;,?/] are ordered with respect to lexicographical ordering. In total degree ordering ()^deg) the monomials 
are ordered with respect to the sum of powers in each monomial. If two monomials have the same sum of powers, 
they are ordered with respect to another ordering, for instance lexicographical. For polynomials in (Q[a;,2/] we have 
x^y ;^deg 4xy since the sum of powers of the left power product is 3 compared to 2 for the right power product. 

The largest power product with respect to the underlying ordering {>-) of a polynomial / is denoted as the leading 
power product, LP(/), the corresponding coefficient as leading coefficient, LC(/). With help of these preparations we 
can define the two essential parts of the Buchberger algorithm, that is Reduction and the S-polynomial. 

Definition. Reduction 

Let f,p (z K[x]. We call f reducible modulo p, if for a power product t of f there exists a power product u with 
LP(p) ■ u — t. Then we say, f reduces to h modulo p, where h — f — c^°''ST^'^^ent(f,t) , , 



In the example (HI I the polynomial /a is reducible modulo /i, since for example the second monomial of /a, that 
is x^y, is a multiple of the LP(/i), and h—fz — l/3/i — x^y — 2/3xy — 1/3?/ + 3x^ — x"^ — 5/3x + 1. 

Reduction of a polynomial modulo a set P C K[x] is accordingly defined if there is a p S P such that / is reducible 
modulo p. Further, we say, a polynomial h is in reduced form or normal form modulo F, in short normf (/i, if 
there is no h' such that h reduces to h' modulo F. A set P C K[x:] is called reduced, if each p € P is in reduced form 
modulo P\{p}. Note that reduction is defined with respect to the underlying ordering of the monomials, since the 
leading power product is defined with respect to the ordering. In general, a normal form is not unique, neither for a 
polynomial nor for a set. 

Now we can present an algorithm, to compute a normal form Q d K[x] of a finite F C /"^^ix]. 

Algorithm. Normal form 

For a given finite set F C K[x] determine a normal form Q C J^[x]. 
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while exists p £ Q 

which is reducible modulo Q\{p} do 

Q Q\{p} 
h := normf (p, Q) 
iih^O then 
Q := Q U {h} 
return Q 



Clearly, the simultaneous zeros of all fiCzF are also simultaneous zeros of all qi £ Q and vice versa. 
Definition. S-polynomial 

For gi,g2 G the S-polynomial of gi and g2 is defined as 

lcm(LP(ffi),LP(g2)) 
spol(ffi,<7,) = 51 

LC(gi)lcm(LP(.gi),LP(g2)) 
LC(52) LP(52) 

where 1cm denotes the least common multiple. 



Clearly, a simultaneous zero of gi and 52 is also a zero of spol((7i, (72)- In the example (HI I we can build the S- 
polynomial for any two polynomials, for instance 

spol(/i, /2) = 4^ A - ^ ^ /2 = X A - 3/2 /2 
x'^y 2 x'^y 

„2„, I c In I o /r)„, I 0^2 



= 2a;^y + b/2xy + 3/2y + Sx^ + 3/2a; - 9/2. 



Finally we define the Groebner basis. 



Definition. Groebner basis 

G C K\x\ is called Groebner Basis, if for all /i, /2 G G normf (spol(/i, A), G) = 0. 

Now we are in a position to present the Buchberger algorithm. 
Algorithm. Buchberger 

For a given finite set F C K[x] determine the Groebner basis G C K[x] with I{F) = I{G). 



G:=F 

B ■■= {{gi,ff2}|5i,52 e G with gi ^ 32} 
while B 7^ do 

choose {51,52} from B 

B :^ i3\{{5i,52}} 

h := spol(5i,52) 

h' := normf (/i, G) 

if ft-V then 

B:=B\j{{g,h'}\g£G} 
G:=GU{/i'} 
return G 
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Note, that since G just follows by adjoining reduced S-polynomials to F both sets generate the same ideal, especially, 
both sets have exactly the same simultaneous zeros. It can be proven, that the Buchberger algorithm terminates. 
The final step is to construct the reduced Groebner basis by applying the normal form algorithm defined above to 
the Groebner basis G. It can be shown that the reduced Groebner basis is unique [207 . If we apply the Buchberger 



algorithm to the set (HI) with subsequent reduction we end up with the reduced Groebner basis (with underlying 



lexicographical ordering) 

gi^y + x'^- i/2x - 3, 
£,2 = -5/2x2 -5/2a;. 

The system of equations 51 = 52 = is equivalent to /i = /2 = /a = 0, but the former allows to directly calculate the 
solutions: Since 172 = is univariate it can be solved immediately and subsequently gi = Q for each partial solution 
inserted. 

Despite the correctness of the Buchberger algorithm tractability of practical examples requires to improve this 
algorithm. In particular, the number of iterations in the algorithm drastically grows with an increasing number of 
polynomials and with higher degrees of the polynomials. In this respect much progress has been made with the 
improvement of this original Buchberger algorithm from 1965; see |207} 12081 1221j . 
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